Pairing computation device, pairing computation method and recording medium storing pairing computation program

ABSTRACT

There is provided a pairing computation device provided with a CPU which computes pairing e(S,Q) with SεG 1 , QεG 2 , χ being a given integer variable, and F being a rational function calculated using Miller&#39;s algorithm with respect to multi-pairing (MMA). An order r, a trace t of the Frobenius endomorphism φ p  are specified preliminarily using the integer variable χ according to an embedding degree k. The CPU performs pairing computation by means of: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a computation unit which computes F; a computation unit which computes a value at a rational point Q(x Q , y Q ) of a straight line passing through given rational points; a computation unit which computes f′  χ,S (Q) using the aforementioned F and the value; and a computation unit which computes the pairing e(S,Q) using aforementioned f′  χ,S (Q) as [F83]. 
     [F83]
 
 e ( S,Q )= f′   χ,S ( Q ) (p     k     −1)/r .

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims priority from PCT/JP2010/057534 filed Apr. 21, 2010, the entire contents of which are incorporated herein by reference, and which in turn claims priority from JP Ser. No. 2009-121398 filed Apr. 21, 2009.

FIELD OF THE INVENTION

The present invention relates to a pairing computation device and a pairing computation method which enable to perform pairing computation at high speed, and a recording medium storing a pairing computation program.

DESCRIPTION OF THE RELATED ART

Conventionally, when an individual user uses various kinds of services provided on a network such as the Internet, there are some cases where authentication procedure using ID and password set preliminarily for each individual user is requested. The authentication procedure like this is for confirming, in general, that an individual user is a normal user and the authentication processing for the authentication procedure is performed by an authentication server.

Recently, with the use of digital signature technology, user-specific digital signature data is added to each individual data per se to make it possible to assure that data used by an individual user is not falsified by a third party or is not leaked to a third party. Accordingly, even high confidential information has been safely handled on a network.

On the other hand, in the digital signature, since an individual user is identified along with authentication processing by an authentication server, each time the authentication processing is performed, a history of each individual user is sequentially accumulated as information on the authentication server. Therefore, since personal information such as what kind of sites did a user get access to, or what kind of services did a user use is accumulated on the authentication server, adequate attention is paid so that any such information should not be leaked from the point of view of personal information protection.

In order to resolve the accumulation of history information about an individual user which arises from employing this digital signature, there has been proposed to employ digital group signature which is extended from the digital signature.

In a case where the digital group signature is employed, an individual user transmits anonymously signature data which certifies only that the individual user belongs to a certain group to an authentication server, and the authentication server authenticates that the individual user belongs to a certain group without identifying individual user from the received signature data. Therefore, the authentication server blocks unfair use by individual user who does not belong to the group on one hand, and authenticates individual user without accumulating history information about each individual user on the other hand.

In the anonymous authentication in the digital group signature, paring computation is employed.

Paring computation is a computation using a function of two inputs and one output. For example, letting S be a rational point over a prime field F_(p), Q be a rational point over a k-th extension field F_(p) ^(k), by inputting two rational points S and Q, an element z over an extension field F*_(p) ^(k) is outputted. And yet, authentication is performed by using a property of bilinearity that when a times of rational point S and b times of rational point Q are inputted, z to the power of ab is outputted. Here, “k” is called an embedding degree, and “F*_(p) ^(k)” is meant to be correctly displayed in mathematical notation as [F31] F*_(p) ^(k) However, due to display restrictions, it is denoted as F*_(p) ^(k).

Generally, for rational points S and Q, points on an elliptic curve are used respectively. The pairing computation of points on an elliptic curve is constituted of a step which computes using Miller's algorithm and a step which performs an exponentiation with respect to a result of the computation.

In the digital group signature, when performing authentication processing of access right of an individual user who belongs to a group, after having performed pairing computation to exclude individual users who have lost access rights, authentication processing is performed by performing pairing computation of the given individual user and hence, it has been possible to respond flexibly to attributes changes due to issuance and expiration of access right for each individual user.

Therefore, in a case of digital group signature where a group is constituted of, for example, 10,000 individual users, if there are 100 individual users who have lost their access rights, it is necessary to perform 100 pairing computations. Since there needs about 0.1 second per one pairing computation performed by present general electronic computer, 100 pairing computations require about 10 seconds. Accordingly, in practical use, the number of individual users is limited and hence, the digital group signature has not been widely used.

Accordingly, in order to enhance usefulness of the digital group signature by increasing computation speed of the pairing computation, there have been proposed techniques which try to realize high-speed computation by, for example, employing Tate pairing computation method which is defined on an elliptic curve as pairing computation to reduce computation load (For example, see JP 2005-316267A.).

SUMMARY OF THE INVENTION

However, speeding up of pairing computation currently proposed is not yet sufficient and further speeding up of the paring computation has been sought.

The inventors, in view of the present situation, have done research and development to speed up pairing computation and have made the present invention.

According to an aspect of the present invention, there is provided a pairing computation device, in which, as described in claim 1, an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), and the pairing computation device comprises a CPU which computes the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, and F be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA),

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified preliminarily according to the embedding degree k using the integer variable χ, and

the CPU includes:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a computation unit which computes F;

a computation unit which computes a value of a straight line passing through given rational points at a rational point Q(x_(Q), y_(Q));

a computation unit which computes f′ _(χ,S)(Q) using said F and said value; and

a computation unit which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F32] e(S, Q)=f′ _(χ,S)(Q)^(p) ^(k) ^(−1)/r).

According to a further aspect of the present invention, there is provided a pairing computation method, in which, as described in claim 2, an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), an electronic computer which includes a CPU computes the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, and F be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA),

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified preliminarily according to the embedding degree k using the integer variable χ, and

the pairing computation method comprises:

a step of inputting the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit;

a step of computing F by causing the CPU of the electronic computer to function as a computation unit;

a step of computing a value of a straight line passing through given rational points at a rational point Q (x_(Q), y_(Q)) by causing the CPU of the electronic computer to function as a computation unit;

a step of computing f′ _(χ,S)(Q) using said F and said value by causing the CPU of the electronic computer to function as a computation unit; and

a step of computing the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F33] e(S, Q)=f′ _(102 ,S)(Q)^((p) ^(k) ^(−1)/r) by causing the CPU of the electronic computer to function as a computation unit.

According to a further aspect of the present invention, there is provided a pairing computation program stored in a recording medium, in which, as described in claim 3,

an elliptic curve is given as y²=x³+ax+b, aε_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), and the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, and F be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA),

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified preliminarily according to the embedding degree k using the integer variable χ, and

the pairing computation program causes the CPU of the electronic computer to function as:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a computation unit which computes F;

a computation unit which computes a value of a straight line passing through given rational points at a rational point Q(x_(Q), y_(Q));

a computation unit which computes f′ _(χ,S)(Q) using said F and said value;

and

a computation unit which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F34] e(S, Q)=f′ _(χ,S)(Q)^((p) ^(k) ^(−1)/r).

According to a further aspect of the present invention, there is provided a pairing computation device, in which, as described in claim 4,

an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ¹²/(F* _(p) ¹²)^(r), and the pairing computation device comprises a CPU which computes the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and pS, Z_(Q) be a set of rational points pQ and Q, and F₂ _(χ,ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA),

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=36χ⁴−36χ³+18χ²−6χ+1, t(χ)=6χ²+1, and a representation of the integer variable χ using p¹⁰ with p as a characteristic is p≡(2χ−1)p¹⁰+2χ(modr(χ)), and the CPU includes:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a first computation unit which computes F_(2χ,ZS)(Z_(Q));

a second computation unit which computes given rational points using 2χS and 2χpS which are calculated when computing said F_(2χ, ZS)(Z_(Q));

a third computation unit which computes a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through the given rational points;

a fourth computation unit which computes f′ _(χ,S)(Q) using said F_(2χ,ZS)(Z_(Q)) and said value; and

a fifth computation unit which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F35] e(S, Q)=f′ _(χ,S)(Q)^((p) ¹² ^(−1)/r).

In particular, the pairing computation device is characterized in that the second computation unit computes respective rational points −S, (2χ−1)S, p¹⁰((2χ−1)S), −pS, (2χ−1)pS, p¹⁰(2χ−1)pS in order using previously obtained results,

the third computation unit respectively computes a value l₁ at a rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2 χ−1)S, −S), a value l₂ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)pS, −pS), a value l₃ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)S), 2χS), and a value l₄ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)pS), 2χpS), and

the fourth computation unit computes f′ _(χ,S)(Q) as

[F36] f′ _(χ,S)(Q)={F _(2χ,Z) _(S) (Z _(Q))·{l ₁ ^(p) ·l ₂}⁻¹}^(p) ₁₀ ·F _(2χ,Z) _(S) (Z _(Q))·l ₃ ^(p) ·l ₄.

According to a further aspect of the present invention, there is provided a pairing computation method, in which, as described in claim 6,

an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ¹²/(F _(p) ¹²)^(r), an electronic computer which includes a CPU computes the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and pS, Z_(Q) be a set of rational points pQ and Q, and F_(2χ,ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA),

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=36χ⁴−36χ³+18χ²−6χ+1 t(χ)=6χ²+1, and a representation of the integer variable χ using p¹⁰ with p as a characteristic is p≡(2χ−1)p ¹⁰+2χ(modr(χ)), and the pairing computation method comprises:

an input step which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit;

a first computation step which computes F_(2χ,ZS)(Z_(Q)) by causing the CPU of the electronic computer to function as a first computation unit;

a second computation step which computes given rational points using 2χS and 2χpS which are calculated when computing said F_(2χ,ZS)(Z_(Q)) by causing the CPU of the electronic computer to function as a second computation unit;

a third computation step which computes a value at a rational point Q (x_(Q), y_(Q)) of a straight line passing through the given rational points by causing the CPU of the electronic computer to function as a third computation unit;

a fourth computation step which computes f′ _(χ,S)(Q) using said F_(2χ,ZS)(Z_(Q)) and said value by causing the CPU of the electronic computer to function as a fourth computation unit; and

a fifth computation step which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F37] e(S, Q)=f′ _(χ,S)(Q)^((p) ¹² ^(−1)/r), by causing the CPU of the electronic computer to function as a fifth computation unit.

In particular, the pairing computation method is characterized in that the second computation step computes respective rational points −S, (2χ−1)S, p¹⁰((2χ−1)S), −pS, (2χ−1)pS, p¹⁰((2χ−1)pS) in order using previously obtained results,

the third computation step respectively computes a value l₁ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2 χ−1)S, −S), a value l₂ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)pS, −pS), a value l₃ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)S), 2χS), and a value l₄ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)pS), 2χpS), and

the fourth computation step computes f′ _(χ,S)(Q) as

[F38] F′ _(χ,S)(Q)={F _(2χ,Z) _(S) (Z _(Q))·{l ₁ ^(p) ·l ₂}⁻¹}^(p) ¹⁰ ·F _(2χ,Z) ^(S) (Z _(Q))·l ₂ ^(p) ·l ₄.

According to a further aspect of the present invention, there is provided a pairing computation program, in which, as described in claim 8, an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ¹²/(F* _(p) ¹²)^(r), the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and pS, Z_(Q) be a set of rational points pQ and Q, and F_(2χ,ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA),

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=36χ⁴−36χ³+18χ²−6χ+1, t(χ)=6χ²+1, and a representation of the integer variable χ using p¹⁰ with p as a characteristic is p≡(2χ−1)p ¹⁰+2χ(modr(χ)), the pairing computation program causes the CPU of the electronic computer to function as:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a first computation unit which computes F_(2χ,ZS)(Z_(Q));

a second computation unit which computes given rational points using 2χS and 2χpS which are calculated when computing said F_(2χ,ZS)(Z_(Q));

a third computation unit which computes a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through the given rational points;

a fourth computation unit which computes f′ _(χ,S)(Q) using said F_(2χ,ZS)(Z_(Q)) and said value; and

a fifth computation unit which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F39] e(S, Q)=f′ _(χ,S)(Q)^((p) ¹² ^(−1)/r).

In particular, the pairing computation program is characterized in that

the pairing computation program causes:

the CPU of the electronic computer which functions as the second computation unit to compute respective rational points −S, (2χ−1)S, p¹⁰((2 χ−1)S), −pS, (2χ−1)pS, p¹⁰((2χ−1)_(PS)) in order using previously obtained results;

the CPU of the electronic computer which functions as the third computation unit to compute respectively a value l₁ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)S, −S), a value l₂ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)pS, −pS), a value l₃ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)S), 2χS), and a value l₄ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)pS, 2χpS); and

the CPU of the electronic computer which functions as the fourth computation unit to compute f′ _(χ,S)(Q) as

[F40] f′ _(χ,S)(Q)={F _(2χ,Z) _(S) (Z _(Q))·{l ₁ ^(p) ·l ₂}⁻¹}^(p) ¹⁰ ·F _(2χ,Z) _(S) (Z _(Q))·l ₃ ^(p) ·l ₄.

According to a further aspect of the present invention, there is provided a pairing computation device, in which, as described in claim 10,

an elliptic curve is given as y²=x³+ax, a∩F_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ⁸/(F* _(p) ⁸)^(r), the pairing computation device comprising a CPU which computes the pairing e (S, Q), by letting SεG, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and p³S, Z_(Q) be a set of rational points p³Q and Q, and F_(3χ,ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA),

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=9χ⁴+12χ³+8χ²+4χ+1, t(χ)=−9χ³−3χ²−2χ, and a representation of the integer variable χ using p² and p³ with p as a characteristic is p³≡p²+3χ+1(modr(χ)) and the CPU includes:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a first computation unit which computes F_(3χ,ZS)(Z_(Q))

a second computation unit which computes respective rational points p²(S), p²(p³S), (3χ+1)S, (3χ+1)p²S in order using previously obtained results;

a third computation unit which respectively computes a value l₅ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χS,S), a value l₆ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3χ+1)S), a value l₇ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χp³S, p²S), and a value l₈ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p² (p₃S), (3χ+1)p³S);

the fourth computation unit which computes f′ _(χ,S)(Q) using a computation result of the first computation unit and a computation result of the third computation unit as

[F41] f′ _(χ,S)(Q)=F _(3χ,Z) _(Z) (Z _(Q)){l ₅ ·l ₆}^(p) ⁵ ·l ₇ ·l ₈; and

a fifth computation unit which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F42] e(S, Q)=f′ _(χ,S)(Q)^((p) ⁸ ^(−1)/r).

According to a further aspect of the present invention, there is provided a pairing computation method, in which, as described in claim 11,

an elliptic curve is given as y²=x³+ax, aεF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ⁸/(F* _(p) ⁸)^(r)), an electronic computer which includes a CPU computes the pairing e (S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and p³S, Z_(Q) be a set of rational points p³Q and Q, and F_(3χ,ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA),

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=9χ⁴+12χ³+8χ²+4χ+1, t(χ)=−9χ³−3χ²−2χ, and a representation of the integer variable χ using p² and p³ with p as a characteristic is p³≡p²+3χ+1(modr(χ)) the pairing computation method comprises:

an input step which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit;

a first computation step which computes F_(3χ,ZS)(Z_(Q)) by causing the CPU of the electronic computer to function as a first computation unit;

a second computation step which computes respective rational points p²(S), p²(p³S), (3χ+1)S, (3χ+1)p³S in order using previously obtained results by causing the CPU of the electronic computer to function as a second computation unit;

a third computation step which respectively computes a value l₅ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χS, S), a value l₆ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3χ+1)S), a value l₇ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χp³S,p³S), and a value l₈ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3χ+1)p³S) by causing the CPU of the electronic computer to function as a third computation unit;

a fourth computation step which computes f′ _(χ,S)(Q) using said F_(3χ,ZS)(Z_(Q)) and said values ₅, l₆, l₇, l₈ as

[F43] f′ ^(χ,S)(Q)=F _(3χ,Z) _(S) (Z _(Q)){l ₅ ·l ₆}^(p) ⁸ ·l ₇ ·l ₈

by causing the CPU of the electronic computer to function as a fourth computation unit; and

a fifth computation step which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F44] e(S, Q)=f′ _(χ,S)(Q)^((p) ⁸ ^(−1)/r) by causing the CPU of the electronic computer to function as a fifth computation unit.

According to a further aspect of the present invention, there is provided a pairing computation program, in which, as described in claim 12,

an elliptic curve is given as y²=x³+ax, aεF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ⁸/(F* _(p) ⁸)^(r), the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and p³S, Z_(Q) be a set of rational points p³Q and Q, and F_(3χ, ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA),

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=9χ⁴+12χ³+8χ²+4χ+1, t(χ)=−9χ³−3χ²−2χ, and a representation of the integer variable χ using p² and p³ with p as a characteristic is p³≡p²+3χ+1(modr(χ)) the pairing computation program causes the CPU of the electronic computer to function as:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a first computation unit which computes F_(3χ, ZS)(Z_(Q));

a second computation unit which computes respective rational points p²(S), p²(p³S), (3χ+1)S, (3χ+1)p³S in order using previously obtained results;

a third computation unit which respectively computes a value l₅ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χS, S), a value l₆ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3χ+1)S), a value l₇ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χp³S, p³S), and a value l₈ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3χ+1)p³S);

a fourth computation unit which computes f′ χ,S(Q) using said F_(3χ,ZS)(Z_(Q)) and said values l₅, l₅, l₇, l₈ as

[F45] f′ _(χ,S)(Q)=F _(3χ,Z) _(S) (Z _(Q)){l ₅ ·l ₆}^(p) ⁵ ·l ₇ ·l ₈; and

a fifth computation unit which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F46] e(S, Q)=f′ _(χ,S)(Q)^(p) ⁸ ^(−1)/r).

According to a further aspect of the present invention, there is provided a pairing computation device, in which, as described in claim 13,

an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(^(SM) _(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), the pairing computation device comprises a CPU which computes the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f be a rational function which is calculated using Miller's algorithm,

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified preliminarily according to the embedding degree k using the integer variable χ, and

the CPU includes:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a computation unit which computes f;

a computation unit which computes a value of a straight line passing through given rational points at a rational point Q(x_(Q), y_(Q));

a computation unit which computes f′ _(χ,S)(Q) using said f and said value; and

a computation unit which computes the pairing e (S, Q) using said f′ _(χ,S)(Q) as

[F47] e(S, Q)=f′ _(χ,S)(Q)^((p) ^(k) ^(−1)/r).

According to a further aspect of the present invention, there is provided a pairing computation method, in which, as described in claim 14,

an elliptic curve is given as y²=x4ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), an electronic computer which includes a CPU computes the pairing e (S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f be a rational function which is calculated using Miller's algorithm,

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified preliminarily according to the embedding degree k using the integer variable χ, and

the pairing computation method comprises:

a step of inputting the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit;

a step of computing f by causing the CPU of the electronic computer to function as a computation unit;

a step of computing a value of a straight line passing through given rational points at a rational point Q(x_(Q), y_(Q)) by causing the CPU of the electronic computer to function as a computation unit;

a step of computing f′ _(χ,S)(Q) using said f and said value by causing the CPU of the electronic computer to function as a computation unit; and

a step of computing the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F48] e(S, Q)=f′ _(χ,S)(Q)^((p) ^(k) ^(−1)/r) by causing the CPU of the electronic computer to function as a computation unit.

According to a further aspect of the present invention, there is provided a pairing computation program in which, as described in claim 15,

an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f be a rational function which is calculated using Miller's algorithm,

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified preliminarily according to the embedding degree k using the integer variable χ, and

the pairing computation program causes the CPU of the electronic computer to function as:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a computation unit which computes f;

a computation unit which computes a value of a straight line passing through given rational points at a rational point Q(x_(Q), y_(Q));

a computation unit which computes f′ _(χ,S)(Q) using said f and said value; and

a computation unit which computes the pairing e (S, Q) using said f′ _(χ,S)(Q) as

[F49] e(S, Q)=f′ _(χ,S)(Q)^((p) ^(k) ^(−1)/r).

According to a further aspect of the present invention, there is provided a pairing computation device, in which, as described in claim 16,

an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ¹²/(F* _(p) ¹²)^(r), the pairing computation device comprises a CPU which computes the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, f_(2χ,S)(Q) be a rational function which is calculated using Miller's algorithm,

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=36χ⁴−36χ³+18χ²−6χ+1, t(χ)=6χ²+1,

a representation of the integer variable χ using p¹⁰ with p as a characteristic is p≡(2χ−1)p¹⁰+2χ(modr(χ)) and the CPU includes:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a first computation unit which computes f_(2χ,S)(Q) and f_(2χ,pS)(Q);

a second computation unit which computes given rational points using 2χS and 2χpS which are calculated when computing said f_(2χ,S)(Q) and f_(2χ,pS)(Q);

a third computation unit which computes a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through the given rational points;

a fourth computation unit which computes f′ _(χ,S)(Q) using said f_(2χ,S)(Q) and said value; and

a fifth computation unit which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F50] e(S, Q)=f′ _(χ,S)(Q)^((p) ¹² ^(−1)/r).

In particular, the pairing computation device is characterized in that

the second computation unit computes respective rational points −S, (2χ−1)S, p¹⁰((2χ−1)S), −pS, (2χ−1)pS, p¹⁰((2χ−1)pS in order using previously obtained results,

the third computation unit respectively computes a value l₁ at a rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2 χ−1)S, −S), a value l₂ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)S ), 2χS), a value l₃ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)pS, −pS), and a value l₄ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1) pS), 2χpS), and

the fourth computation unit computes f′ _(χ,S)(Q) using the values of the rational point Q (x_(Q), y_(Q)) l₁, l₂, l₃, l₄ as

[F51] f′ _(χ,S)(Q)=({f _(2χ,S)(Q)·l ₁ ⁻¹}^(p) ¹⁰ ·f _(2χ,S)(Q)·l ₂)^(p) ·{f _(2χ,pS)(Q)·l ₃ ⁻¹}^(p) ¹⁰ ·f _(2χ,pS)(Q)·l ₄.

According to a further aspect of the present invention, there is provided a pairing computation method, in which, as described in claim 18,

an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ¹²/(F* _(p) ¹²)^(r), an electronic computer which includes a CPU computes the pairing e (S, Q), by letting SεG_(I), QεG₂, χ be a given integer variable, and f_(2χ,S)(Q) be a rational function which is calculated using Miller's algorithm,

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=36χ⁴−36χ³+18χ²−6χ+1, t(χ)=6χ²+1,

a representation of the integer variable χ using p¹⁰ with p as a characteristic is p≡(2χ−1)p¹⁰+2χ(modr(χ)), and

the pairing computation method comprises:

an input step which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit;

a first computation step which computes f_(2χ,S)(Q) and f_(2χ,pS)(Q) by causing the CPU of the electronic computer to function as a first computation unit;

a second computation step which computes given rational points using 2χS and 2χpS which are calculated when computing said f_(2χ,S)(Q) and f_(2χ,pS)(Q) by causing the CPU of the electronic computer to function as a second computation unit;

a third computation step which computes a value at a rational point Q (x_(Q), y_(Q)) of a straight line passing through the given rational points by causing the CPU of the electronic computer to function as a third computation unit;

a fourth computation step which computes f′ _(χ,S)(Q) using said f_(2χ,S)(Q) and said value by causing the CPU of the electronic computer to function as a fourth computation unit; and

a fifth computation step which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F52] e(S, Q)=f′ _(χ,S)(Q)^((p) ¹² ^(−1)/r), by causing the CPU of the electronic computer to function as a fifth computation unit.

In particular, the pairing computation method is characterized in that

the second computation step computes respective rational points −S, (2χ−1)S, p¹⁰((2χ−1)S), −pS, (2χ−1)pS, p¹⁰((2χ−1)pS in order using previously obtained results,

the third computation step respectively computes a value l₁ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2 χ−1)S, −S), a value l₂ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)S), 2χS), a value l₃ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1) pS, −pS), and a value l₄ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)pS), 2χpS), and

the fourth computation step computes f′ _(χ,S)(Q) using the values of the rational point Q (x_(Q), y_(Q)) l₁, l₂, l₃, l₄ as f′ _(χ,S)(Q)=({f _(2χ,S)(Q)·l ₁ ⁻¹}^(p) ¹⁰ ·f _(2χ,S)(Q)·l ₂)^(p) ·{f _(2χ,pS)(Q)·l ₃ ⁻¹}^(p) ¹⁰ ·f _(2χ,pS)(Q)·l ₄.

According to a further aspect of the present invention, there is provided a pairing computation program, in which, as described in claim 20,

an elliptic curve is given as y²=x³b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ¹²/(F* _(p) ¹²)^(r), and the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f_(2χ,S)(Q) ^(and f) _(2χ,pS)(Q) be a rational function which is calculated using Miller's algorithm,

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=36χ⁴−36χ³+18χ²−6χ+1, t(χ)=6χ²+1,

a representation of the integer variable χ using p¹⁰ with p as a characteristic is p≡(2χ−1)p¹⁰+2χ(modr(χ)), and

the pairing computation program causes the CPU of the electronic computer to function as:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a first computation unit which computes f_(2χ, S)(Q) and f_(2χ,ps)(Q);

a second computation unit which computes given rational points using 2χS and 2χpS which are calculated when computing said f_(2χ,S)(Q) and f_(2χ,pS)(Q);

a third computation unit which computes a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through the given rational points;

a fourth computation unit which computes f′ _(χ,S)(Q) using said f_(2χ, S)(Q) and said value; and

a fifth computation unit which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F54] e(S, Q)=f′ _(χ,S)(Q)^((p) ¹² ^(−1)/r).

In particular, the pairing computation program is characterized in that

the pairing computation program causes:

the CPU of the electronic computer which functions as the second computation unit to compute respective rational points —S, (2χ−1)S, p¹⁰((2 χ−1)S), −pS, (2χ−1)pS, p¹⁰((2χ−1)pS) in order using previously obtained results;

the CPU of the electronic computer which functions as the third computation unit to respectively compute a value l₁ at a rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)S, −S), a value l₂ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)S), 2χS), a value l₃ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)pS, −pS), and a value l₄ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)pS), 2χpS); and

the CPU of the electronic computer which functions as the fourth computation unit to compute f′ _(χ,S)(Q) using the values at the rational point Q (x_(Q), y_(Q)) l₁, l₂, l₃, l₄ as

[F55] f′ _(χ,S)(Q)=({f _(2χ,S)(Q)·l ₁ ⁻¹}^(p) ¹⁰ ·f _(2χ,S)(Q)·l ₂)^(p) ·{f _(2χ,pS)(Q)·l ₃ ⁻¹}^(p) ¹⁰ ·f _(2χ,pS)(Q)·l ₄.

According to a further aspect of the present invention, there is provided a pairing computation device, in which, as described in claim 22,

an elliptic curve is given as y²=x³+ax, aεF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ⁸/(F* _(p) ⁸)^(r), and the pairing computation device comprises a CPU which computes the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f_(3χ,S)(Q) be a rational function which is calculated using Miller's algorithm,

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=9χ⁴+12χ³+8χ²+4χ+1, t(χ)=−9χ³−3χ²−2χ,

a representation of the integer variable χ using p² and p³ with p as a characteristic is p³≡p²+3χ+1(modr(χ)) and the CPU includes:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a first computation unit which computes f_(3χ,S)(Q) and f_(3χ,p^3S)(Q);

a second computation unit which computes respective rational points p²(S), (3χ+1)S, p²(p³S), (3χ+1)p³S in order using previously obtained results;

a third computation unit respectively computes a value l₅ at a rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (3χS, S), a value l₆ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3χ+1)S), a value l₇ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (3χp³S, p³S), and a value l₈ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p² (p³S), (3χ+1)p³S); and

the fourth computation unit which computes f′ _(χ,S)(Q) using the computation result of the first computation unit and the computation result of the third computation unit as

[F56 ] f′ _(χ,S)(Q)=(f _(3χ,S)(Q)·l ₅ ·l ₆)^(p) ⁵ f _(3χ,p) 3 _(S)(Q)·l ₇ ·l ₈; and

a fifth computation unit which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F57] e(S, Q)=f′ _(χ,S)(Q)^((p) ⁹ ^(−1)/r).

According to a further aspect of the present invention, there is provided a pairing computation method, in which, as described in claim 23,

an elliptic curve is given as y²=x³+ax, aÅF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ⁸/(F* _(p) ⁹)^(r), an electronic computer which includes a CPU computes the pairing e(S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f_(3χ, S)(Q) be a rational function which is calculated using Miller's algorithm,

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=9χ⁴+12χ³+8χ²+4χ+1, t(χ)=−9χ³−3χ²−2χ,

a representation of the integer variable χ using p² and p³ with p as a characteristic is p³≡p²+3χ+1(modr(χ)) and

the pairing computation method comprises:

an input step which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit;

a first computation step which computes f_(3χ,S)(Q) and f_(3χ,p^3S)(Q) by causing the CPU of the electronic computer to function as a first computation unit;

a second computation step which computes respective rational points p²(S), (3χ+1)S, p²(p³S), (3χ+1)p³S in order using previously obtained results by causing the CPU of the electronic computer to function as a second computation unit;

a third computation step which respectively computes a value l₅ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χS, S), a value l₆ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3χ+1)S), a value l₇ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χp³S, p³S), and a value l₈ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3χ+1)p³S) by causing the CPU of the electronic computer to function as a third computation unit; and

a fourth computation step which computes f′ _(χ,S)(Q) using the values l₅, l₆, l₇, l₈ as

[F58] f′ _(χ,S)(Q)=(f _(3χ,S)(Q)·l ₅ ·l ₆)^(pds 3) f _(3χ,p3) _(S)(Q)·l ₇ ·l ₈;

by causing the CPU of the electronic computer to function as a fourth computation unit; and

a fifth computation step which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F59] e(S, Q)=f _(χ,S)(Q)^((p) ⁸ ^(−1)/r) by causing the CPU of the electronic computer to function as a fifth computation unit.

According to a further aspect of the present invention, there is provided a pairing computation program, in which, as described in claim 24,

an elliptic curve is given as y²=x³+ax, aεF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ⁸/(F* _(p) ⁸)^(r), the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e (S, Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f_(3χ,S)(Q) be a rational function which is calculated using Miller's algorithm,

the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=9χ⁴+12χ³+8χ²+4χ+1, t(χ)=−9χ³−3χ²−2χ,

a representation of the integer variable χ using p² and p³ with p as a characteristic is p³≡p²+3χ+1(modr(χ)) and

the pairing computation program causes the CPU of the electronic computer to function as:

an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers;

a first computation unit which computes f_(3χ,S)(Q) and f_(3χ,p^S)(Q);

a second computation unit which computes respective rational points p²(S), (3χ+1)S, p²(p³S), (3χ+1)p′S in order using previously obtained results;

a third computation unit which respectively computes a value l₅ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χS, S), a value l₆ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3χ+1)S), a value l₇ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χp³S, p³S), and a value l₈ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3χ+1)p³S); and

a fourth computation unit which computes f′ _(χ,S)(Q) using the values l₅, l₆, l₇, l₈ as

[F60] f′ _(χ,S)(Q)=(f _(2χ,S)(Q)·l ₅ ·l ₆)^(p) ⁸ f _(3χ,p) ₈ _(S)(Q)·l ₇ ·l ₈; and

a fifth computation unit which computes the pairing e(S, Q) using said f′ _(χ,S)(Q) as

[F61] e(S, Q)=f _(χ,S)(Q)^((p) ⁸ ^(−1)/r).

According to the present invention, it is possible to calculate a rational function at high speed by making the rational function a function of an integer variable χ, thus enabling high speed pairing computation.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic diagram of a pairing computation device according to an embodiment of the present invention;

FIG. 2 is a flowchart of a Twisted Xate pairing computation program according to the embodiment of the present invention;

FIG. 3 is a flowchart for computing a rational function f_(2χ,S)(Q);

FIG. 4 is a flowchart of the Twisted Xate pairing computation program according to the embodiment of the present invention;

FIG. 5 is a flowchart of the Twisted Xate multi-pairing computation program according to the embodiment of the present invention;

FIG. 6 is a flowchart for computing a rational function F_(2χ,ZS)(Z_(Q));

FIG. 7 is a flowchart of a Twisted Xate multi-pairing computation program according to the embodiment of the present invention; and

FIG. 8 is a calculation flow in a case where a thread calculation is applied to the Twisted Xate pairing according to the embodiment of the present invention.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

In a pairing computation device, a pairing computation method and a pairing computation program according to the present invention, out of a first step which computes a rational function using Miller's algorithm in pairing computation and a second step which performs exponentiation of a result of the computation, the first step computes the rational function using an integer variable χ, thereby speeding up the pairing computation.

That is, in a conventional pairing computation, in which

an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), a pairing α(S, Q) is computed, using SεG₁, QεG₂, a given integer variable χ, and a rational function f_((t−1)^2,S)(Q) which is calculated using Miller's algorithm, as [F62] α(S, Q)=f _((t−1))Z,_(S)(Q)^((p) ^(k) ^(−1)/r). This is known as Twisted Ate pairing.

On the other hand, the inventors have found out a new pairing which enables fast pairing computation as a non-degenerate bilinear map by using an integer variable χ on an elliptic curve as below e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r). This pairing is hereinafter referred to as “Twisted Xate pairing” . Further, it is possible to perform further faster pairing computation by applying multi-pairing technique which calculates a plurality of pairing products. This pairing is hereinafter referred to as “Twisted Xate multi-pairing” That is, the pairing computation device, the pairing computation method, and the pairing computation program all according to the present invention make it possible to perform fast pairing computation by using Twisted Xate pairing and Twisted Xate multi-pairing.

In particular, elliptic curves have been known as pairing friendly curves according to embedding degrees respectively. For example, when an embedding degree is 12, there has been known that an order r and a trace t of the Frobenius endomorphism φ_(p) are expressed using a integer variable χ as r(χ)=36χ⁴−36χ³+18χ²−6χ+1, t(χ)=6χ²+1.

Further, when the embedding degree is 8, there has been known that r (χ) and t (χ) are expressed as r(χ)=9χ⁴+12χ³+8χ²+4χ+1, t(χ)=−9χ³−3χ²−2χ.

Hereinafter, a case where an embedding degree is 12 is explained as an example. Accordingly, as mentioned above, r(χ)=36χ⁴−36χ³+18χ²−6χ+1   (Equation 1) t(χ)=6χ²+1   (Equation 2) From (Equation 2), 6χ²≡t−1≡p(mod r)   (Equation 3) Here, it is used that a characteristic p is expressed as p=r+t−1   (Equation 4) Accordingly, (Equation 1) is transformed as p²−6χp+3p−6χ+1≡0(mod r) (−6χ+3)p≡−p²+6χ−1(mod r)   (Equation 5) By squaring both sides of (Equation 5), the following equation is obtained. (6χ−3)²p²≡(p²−6χ+1)²(mod r) 36χ²p²−36χp²+9p²≡p⁴−12χp²+2p²+36χ²−12χ+1(mod r)   (Equation 6) Here, from already known relation equation p⁴+1≡p²(mod r), 36χ²p²−36χp³+9p²≡−12χp²+3p²+36χ²−12χ(mod r), 36χ²(p²−1)≡(24χ−6)p²−12χ(mod r), 6χ²(p²−1)≡(4χ−1)p²−2χ(mod r)   (Equation 7) are obtained and by multiplying both sides of (Equation 7) by (p²−1)⁻¹ 6χ²≡−(4χ−1)p⁴+2χp²   (Equation 8) is obtained. Here, (p²−1)⁻¹ can be obtained based on relation equation of p⁴−p²+1≡0(mod r) and gcd(p⁴−p²+1, p²−1)=1, as follows. p⁴−p²+1≡0(mod r) −p²(p²−1)≡1(mod r) (p²−1)⁻¹≡−p²(mod r)   (Equation 9) From the relation equation p²≡p⁴+1(mod r),

$\begin{matrix} \begin{matrix} {{6x^{2}} \equiv {{{- \left( {{4x} - 1} \right)}p^{4}} + {2{x\left( {p^{4} + 1} \right)}}}} \\ {\equiv {{{- \left( {{2x} - 1} \right)}p^{4}} + {2x\;\left( {{mod}\; r} \right)}}} \end{matrix} & \left( {{Equation}\mspace{14mu} 10} \right) \end{matrix}$ is obtained. Based on (Equation 3) and p⁶≡−1(mod r), p≡(2χ−1)p¹⁰+2χ  (Equation 11) is obtained.

Next, a rational function f_((t−1)^2 mod r,S)(·) of Twisted Ate pairing is considered. In particular, from (Equation 3),

[F63] f _(p) _(Z) _(mod r,S) =f _((t−1)) _(Z) _(mod r,S)   (Equation 12) Here, SεG₁ and with respect to QεG₂, [F64] f _(p) _(Z) _(mod r,S)(Q)^((p) ¹² ^(−1)/r) =f _((t−1) _(Z) _(,S)(Q)^((p) ¹² ^(−1)/r)=α(S, Q)   (Equation 13)

Further, a rational function has relationships as follows. f _(a+b,S) =f _(a,S) ·f _(b,S) ·g _(aS,bS),   (Equation 14a) f _(a−b,S) =f _(a,S) ·{f _(b,S) ·g _((a−b)S,bS)}⁻¹   (Equation 14b) f _(ab,S) =f ^(a) _(b,S) ·f _(a,bS) =f ^(b) _(a,S) ·f _(b,aS)   (Equation 14c) Here, g_(aS,bS)=_(aS,bS)/v_(aS+bS), and l_(aS,bS) is a value of a straight line passing through two rational points aS and bS, and v_(aS+bS) is a value of a vertical line at a rational point aS+bS. In a case of the embedding degree being an even number, calculation of v_(aS+bS) can be omitted. Based on (Equation 14c), (Equation 12) becomes as follows. [F65] f _(p) _(Z) _(S) =f _(p,S) ^(p) ·f _(p,pS)   (Equation 15) And using (Equation 11), it becomes [F66] f _(p,S) ^(p) ·f _(p,pS) =f _((2χ−1)p) ₁₀ _(+2χ,S) ^(p) ·f _((2χ−1)p) ₁₀ _(+1χ,S)   (Equation 16) Further, based on (Equation 14a) and (Equation 14b), it becomes [F67] f _(p,S) ^(p) ·f _(p,pS)=({f _(2χ,S) ·{g _((2χ−1)S,−S)}⁻¹}^(p) ¹⁰ ·f _(p) ₁₀ _(,(2χ−1)S) ·f _(2χ,S) ·g _((2χ−1)p) ₁₀ _(,2χS))^(p) ·{f _(2χ,pS) ·{g _((2χ−1)p) ₁₀ _(S,2χS)}⁻¹}^(p 10) ·f _(p) ₁₀ _(,(2χ−1)pS) ·f _(2χ.pS) ·g _((2χ−1)p) ₁₀ _(pS,2χpS)   (Equation 17) In addition, since in (Equation 17), [F68] f _(p) ₁₀ _(,(2χ−1)) ^(p) ·f _(p) ₁₀ _(,(2χ−1)pS)   (Equation 18) has a property of bilinearity, (Equation 17) can be transformed as follows. [F69] f _(p,S) ^(p) ·f _(p,pS) ·{f _(p) ₁₉ _(,(2χ−1)S) ^(p) ·f _(p) ₁₀ _(,(2χ−1)pS)}⁻¹=({f _(2χS) ·{g _((2χ−1)S,−S)}⁻¹}^(p) ¹⁰ ·f _(2χ,S) ·g _((2χ−1)p) ₁₀ _(S,2χS))^(p) ·{f _(2χ,pS) ·{g _((2χ−1)pS,−pS)}⁻¹}^(p) ¹⁰ ·f _(2χ,pS) ·g _((2χ−1) p) ¹⁰ _(pS,2χpS)   (Equation 19) Here, that the left hand side of (Equation 19) has a property of bilinearity brought the inventors to the fact that the right hand side of (Equation 19) also has a property of bilinearity. Accordingly, the inventors have found out Twisted Xate pairing which, using the right hand side of (Equation 19) as a new rational function f′ _(χ,S)(·), computes pairing e(S, Q) as [F70] e(S, Q)=f′ _(χ,S)(Q)^((p) ¹² ^(−1)/r)   (Equation 20). In this Twisted Xate pairing, it is possible to obtain f′ _(χ,S)(Q) at high speed by calculating f_(2χ,S)(Q) and f_(2χ,pS)(Q) using an integer variable χ which has a smaller size than (t−1)² mod r in a computation using Miller's algorithm. However, since two rational points used in these two computations using Miller's algorithm are different respectively, usually (Equation 19) is calculated as follows.

Given F₁, F₂ as

[F71] F ₁ {f _(2χ,S) ·{g _((2χ−1)S,−S)}⁻¹}^(p) ¹⁰ ·f _(2χS) ·g _((2χ−1)p) ₁₀ _(S,2χS)   (Equation 21) [F72] F ₂ ={g _(2χ,pS) ·{g _((2χ−1)pS,−pS)}⁻¹}^(p) ¹⁰ ·f _(2χ,pS) ·g _((2χ−1)p) ¹⁰ _(pS,2χps)   (Equation 22) (Equation 19) can be expressed as f′ _(χ,S)=F₁ ^(p)·F₂ using F₁ and F₂. Since rational point pS which is necessary for calculating F₂ is obtained in the course of calculating F₁, F₁ is calculated first and next, F₂ is calculated. In this case, a repetition count of Miller's algorithm becomes 2[ log_(2χ)]≈(½)[ log₂r].

A basic idea of the present invention is to make F₁ and F₂ mutually independent calculations by calculating pS preliminarily. Due to this, Twisted Xate pairing according to the present invention can realize parallel processing. To be more specific, in an environment where a multi-pairing which enables parallel processing of a plurality of Miller's algorithm calculations and a CPU having a plurality of cores can be used, thread calculation can be applicable.

<Application of Thread Calculation>

pS can be calculated at high speed by means of a scalar multiplication to an extent of (¼) [ log₂r] from (Equation 11).

By preliminarily calculating pS, F₁ and F₂ become mutually independent calculations. Accordingly, in an environment where a CPU having a plurality of cores can be used, parallel processing can be realized by assigning calculations of F₁ and F₂ to respective different cores. FIG. 8 shows a calculation flow in a case where a thread calculation is applied. In this case, a repetition count of Miller's algorithm substantially becomes (¼) [ log₂r] from (Equation 21) and (Equation 22). <Application of Multi-Pairing>

Further, the right hand side of (Equation 19) can be transformed into {f _(2χ,S) ^(p) ·f _(2χ,pS) ·{g _((2χ−1)S,−S) ^(p) ·g _((2χ−1)pS,−pS)}⁻¹}^(p) ¹⁰ ·f _(2χ,S) ^(p) ·f _(2χ,pS) ·g _((2χ−1)p) ₁₀ _(S,2χS) ^(p) ·g _((2χ−1)pS,2χpS)   (Equation 23) From (Equation 23), it is realized that the main calculation of Twisted Xate pairing is A=f^(p) _(2χ,S)·f_(2χ,pS). In general, to calculate A, f_(2χ,S) is calculated first using Miller's algorithm a repetition count of which is given as 2χ, and then, f_(2χ,pS) is calculated. And, A is obtained by taking a product of the two. Here, f^(p) _(2χ,S) can be obtained at an extremely small calculation cost by means of the Frobenius endomorphism. Therefore, since Miller's algorithm a repetition count of which is given as (¼) [ log₂r] is performed twice, the calculation cost of processing of Twisted Xate pairing by means of Miller's algorithm is given as (½) [ log₂r]. The present invention is characterized in that A is obtained by means of Miller's algorithm with respect to multi-pairing (MMA), although in general, A is obtained by independently calculating Miller's algorithm. Using the next relational expression, [F74] f _(r,S)(Q)^(p) =f _(r,S)(pQ), A becomes as below. [F75] A=f _(2χ,S)(pQ)·f _(23χ,pS)(Q). By giving A in this way, it can be thought that the calculation of A is a multi-pairing which calculates a product of two pairings. Accordingly, by applying MMA, (Equation 23) can be calculated as follows. [F76] f′ _(χ,S)(Q)={F _(2χ, Z) _(S) (Z _(Q))·{g_((2χ−1)S,−S) ^(p) ·g _((2χ−1)pS,−pS)}⁻¹}^(p) ¹⁰ ·F _(2χ,Z) _(S) (Z _(Q))·g _((2χ−1)p) ₁₀ _(S,2χS) ^(p) ·g _((2χ−1)pS,2χpS)   (Equation 26) Here, F_(2χ,ZS) is obtained by giving 2χ to a repetition count of MMA algorithm and by giving Z_(S) and Z_(Q) the following equations respectively. Z_(S)={S, pS}  (Equation 27) Z_(Q)={pQ, Q}  (Equation 28) By calculating in this way, it becomes possible to reduce the number of square operations over F_(p) ¹² in Miller's algorithm and the number of inverse element calculations over F_(p) required for elliptic addition and elliptic doubling to ½ as compared with a case where the calculations by means of Miller's algorithm are performed independently. Here, from (Equation 27) and (Equation 28), pS and pQ are required to be preliminarily calculated for applying MMA to the calculation of A. Since pS can be calculated at high speed by means of scalar multiplication to an extent of (¼) [ log₂r] from (Equation 11) and pQ is obtained by means of the Frobenius endomorphism, a calculation cost required for applying the present invention is extremely smaller than a calculation cost which can be reduced by applying MMA.

Heretofore, an explanation is made in a case where an embedding degree is 12. An explanation in a case where the above mentioned embedding degree is 8 is basically the same and hence, detailed explanation is omitted.

Hereinafter, an embodiment of the Twisted Xate pairing in a case where the embedding degree is 12 is explained in detail. In addition, in the present embodiment, digital group signature is assumed and an authentication server constituted of required electronic computer is assumed to be a pairing computation device. However, a pairing computation device is not limited to a case where the pairing computation device is constituted of an authentication server. Provided a device includes at least an arithmetic means such as a CPU and can perform pairing computation, any device may be used as the pairing computation device.

As shown in FIG. 1, an electronic computer 10 which constitutes an authentication server includes: a CPU 11; a storage unit 12 such as a hard disk which stores various kinds of programs such as a pairing computation program, data which is used by the pairing computation program and the like; a memory unit 13 which expands the pairing computation program and makes it executable, temporarily stores data generated in the course of execution of the pairing computation program and is constituted of such as a RAM. In FIG. 1, numeral 14 is a bus.

Further, the electronic computer 10 is connected to a telecommunication line 20 such as the Internet and is configured to be able to receive signature data of digital group signature transmitted from a client device 30 which is connected the telecommunication line 20. In FIG. 1, numeral 15 is an input output control part of the electronic computer 10.

The electronic computer 10, when signature data of the digital group signature is transmitted from the client device 30, stores temporarily the transmitted signature data in the memory unit 13, starts a pairing computation program and performs pairing computation.

That is, the electronic computer 10, by means of the started pairing computation program, performs pairing computation based on a flowchart shown in FIG. 2 and realizes the digital group signature. Here, an authentication processing in the digital group signature is not explained in detail and only the pairing computation as a subroutine processing in the authentication processing is explained in detail.

The electronic computer 10, causing the CPU 11 to function as an input unit by means of the started pairing computation program, inputs data of an integer variable χ and data of a rational point S both of which are stored preliminarily in the memory unit 13, and data of a rational point Q temporarily stored in the memory unit 13 as the signature data into predetermined registers provided inside of the CPU 11 respectively. (Step S1).

Next, the electronic computer 10, causing the CPU 11 to function as a first computation unit by means of the pairing computation program, computes a rational function f_(2χ,S)(Q) by means of Miller's algorithm and stores a result of the computation in the memory unit 13. (Step S2).

The electronic computer 10 performs the computation of the rational function f_(2χ,S)(Q) specifically as shown in FIG. 3. In particular, in step S2, the electronic computer 10 computes 2χS as well as the rational function f_(2 χ,S)(Q) and stores results of the computation in the memory unit 13.

That is, firstly, the electronic computer 10, as initial setting, sets f to be 1 and T to be S and also sets the number of bits in a case of binary representation of the integer variable 2χ to be i (Step S21).

Next, the electronic computer 10 performs predetermined computation of the rational function f_(2χ,S)(Q) part in step S22, performs predetermined computation of 2χ S part in step S23, determines whether u_(i) which is a value of i-th bit counted from the least degree of the integer variable 2χ is “1” or “0” (Step S24) and, when u_(i)=1, further performs predetermined computation of f_(2χ,S)(Q) part in step S25 and performs predetermined computation of 2χS part in step S26.

Next, the electronic computer 10, when (Step S27), decrements i (Step S28), returns to step S22, by repeating the processing until i=1 in step S27 computes the rational function f_(2χ,S)(Q) and 2χS, and stores results of the computation in the memory unit 13 respectively.

Next, the electronic computer 10, causing the CPU 11 to function as a second computation unit by means of the pairing computation program, reads a value of 2χS stored in the memory unit 13, performs computation of rational points −S, (2χ−1)S, p¹⁰((2χ−1)S) respectively and stores values of the respective rational points in the memory unit 13 respectively. (Step S3).

Next, the electronic computer 10, causing the CPU 11 to function as a third computation unit by means of the pairing computation program, reads values of the respective rational points obtained in step S3 from the memory unit 13 respectively, computes l₁ which is a value at a rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)S, −S) and l₂ which is a value at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)S), 2χS) respectively, and stores values of l₁ and l₂ in the memory unit 13 respectively. (Step S4).

To be more specific, the electronic computer 10 computes l=l_(A,B)(Q) as λ_(A,B)←(y_(B)−y_(A))/(x_(B)−x_(A)), l_(A,B)(Q)←(x_(Q)−x_(B))λ_(A,B)−(y_(Q)−y_(B)).

That is, the CPU 11 computes a gradient λ_(A,B) using xy coordinates of points A and B and computes l_(A,B)(Q) using a result of the computation and xy coordinates of points Q and B, and stores a result of the computation in the memory unit 13.

Next, the electronic computer 10, causing the CPU 11 to function as the second computation unit by means of the pairing computation program, computes rational point pS (Step S5). To be more specific, the electronic computer 10 reads values of rational points obtained in step S2 and step S3 from the memory unit 13 and computes using these values as below and stores the value of pS in the memory unit 13. pS=p ¹⁰((2χ−1)S)+2χS

Next, the electronic computer 10, causing the CPU 11 to function as the first computation unit by means of the pairing computation program, computes the rational function f_(2χ,pS)(Q) by means of Miller's algorithm and stores a result of the computation in the memory unit 13. (Step S6).

The electronic computer 10 performs a computation of the rational function f_(2χ,S)(Q) specifically as shown in FIG. 3. In particular, in step S6, the electronic computer 10 performs computation of 2χpS as well as computation of the rational function f_(2χ,S)(Q) and stores a result of the computation in the memory unit 13.

Next, the electronic computer 10, causing the CPU 11 to function as the second computation unit by means of the pairing computation program, reads a value of 2χpS stored in the memory unit 13, performs computation of rational points −pS, (2χ−1)pS, p¹⁰((2χ−1)pS) respectively and stores values of the respective rational points in the memory unit 13 respectively. (Step S7).

Next, the electronic computer 10, causing the CPU 11 to function as the third computation unit by means of the pairing computation program, reads values of the respective rational points obtained in step S7 from the memory unit 13 respectively, computes l₃ which is a value at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)pS,−pS) and l₄ which is a value at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)pS), 2χpS) respectively, and stores values of l₃ and l₄ in the memory unit 13 respectively. (Step S8).

Next, the electronic computer 10, causing the CPU 11 to function as a fourth computation unit by means of the pairing computation program, reads a result of the computation in the first computation unit (values of f_(2χ,S)(Q) and f_(2χ,S)(Q)) and a result of the computation in the third computation unit (values of l₁, l₂, and l₃, l₄) from the memory unit 13, computes f′ _(χ,S)(Q) expressed as below, and stores a result of the computation in the memory unit 13 (Step S9).

[F77] f′ _(χ,S)(Q)=({f _(2χ,S) ·l ₁ ⁻¹}^(p) ¹⁰ ·f _(2χ,S) ·l ₂)^(p) ·{f _(2χ,pS) ·l ₃ ⁻¹}^(p) ¹⁰ ·f _(χ,pS) ·l ₄   (Equation 29)

To be more specific, the electronic computer 10 computes as follows.

input: SεG₁, QεG₂, 2 χ, p output: f=f’ _(χ,) _(S)(Q) procedure:  1. A←f_(2χ,) _(S) //T₁←2 χ S is obtained  2. T₂← −S  3. T₃←T₁ +T₂  4. B←g_(T3, T2)  5. B←A·B⁻¹  6. B←B

p¹⁰  7. f←B·A  8. T₄←p¹⁰T₃  9. B←g_(T4, T1) 10.  T₄←T₄ +T₁ //T₄←pS is obtained 11.  f←f·B 12.  f←f^(p) 13.  A←f_(2χ,) _(T4) //T₁←2 χ pS is obtained 14.  T₂← −T₄ 15.  T₃←T₁ +T₂ 16.  B←g_(T3, T2) 17.  C←A·B⁻¹ 18.  C←C

p¹⁰ 19.  C←C·A 20.  T₄←p¹⁰T₃ 21.  B←g_(T4, T1) 22. C←C·B 23.  f←f·C 24. Return f.

That is, the CPU 11 performs respective computations in the procedure described above using the memory unit 13, obtains f′ _(χ,S)(Q), and stores the value in the memory unit 13. To be more specific, in step 1, the CPU 11 obtains f_(2χ,S) by means of Miller's algorithm and assigns the result to A and also assigns 2χS obtained simultaneously to T₁. In step 2 and 3, the CPU 11 obtains rational points −S and (2χ−1)S, and assigns the results to T₂ and T₃ respectively. In step 4, as mentioned above, since in a case where an embedding degree is an even number, G_(A,B) becomes l_(A,B), the CPU 11 obtains l₁ which is a value at rational point Q of a straight line passing through rational points ((2_(χ−1))S, −S), and assigns the result to B. In step 5 and 6, the CPU 11 obtains a { } part in ( ) of (Equation 29) and assign the result to B. In step 7, the CPU 11 obtains a product of B and A, and assigns the result to f. In step 8, the CPU 11 obtains rational point p¹⁰((2χ−1)S) and assigns the result to T₄. In step 9, the CPU 11 obtains l₂ which is a value at rational point Q of a straight line passing through rational points (p¹⁰((2χ−1)S), 2χS) and assigns the result to B. In step 10, the CPU 11 obtains pS using relationship in (Equation 11) and assigns the result to T₄. In step 11 and 12, the CPU 11 obtains a ( ) part of (Equation 29) and assigns the result to f. In step 13, the CPU 11 obtains f_(2χ,pS) by means of Miller's algorithm and assigns the result to A, and assigns 2χpS obtained simultaneously to T₁. In step 14 and 15, the CPU 11 obtains rational points −pS and (2χ−1)pS. In step 16, the CPU 11 obtains l₃ which is a value at rational point Q of a straight line passing through rational points ((2χ−1)pS, −pS) and assigns the results to B. In step 17 and 18, the CPU 11 obtains a { } part accompanying ( ) in (Equation 29) described above. In step 19, the CPU 11 obtains a product of C and A, and assigns the result to C. In step 20, the CPU 11 obtains rational point p¹⁰((2χ−1)S) and assigns the result to T₄. In step 21, the CPU 11 obtains l₄ which is a value at rational point Q of a straight line passing through rational points (p¹⁰((2χ−1)S), 2χpS) and assigns the result to B. In step 22 and 23, the CPU 11 obtains f′ _(χ,S)(Q) and assigns the result to f.

Next, the electronic computer 10, causing the CPU 11 to function as a fifth computation unit by means of the pairing computation program, reads a value of f′ _(χ,S)(Q) obtained in step S9 from the memory unit 13, performs an exponentiation of the final exponentiation in pairing e(S, Q), and stores a result of the computation in the memory unit 13 (Step S10).

To be more specific, the electronic computer 10 computes as follows and stores f′, that is, a value of e(S, Q) in the memory unit 13.

1. f’ ←f’ 

p⁶ ·f’ ⁻¹ 2. f’ ←f’ 

p² ·f’ 3. a←(f’ ⁶ ) ^(χ) · (f’ ⁵) 

p⁶ 4. b←a^(p) 5. b←a·b 6. compute f’ 

p, f’ 

p², and f’ 

p³ 7. c←b· (f’ 

p)² ·f’ 

p² 8. f’ ←f’ 

p³ · (c⁶) 

 χ² ·c·b· (f’ 

p ·f’ )⁹ ·a·f’ 

 4 9. Return f’

The electronic computer 10 which constitutes an authentication server performs an authentication processing using a result of the pairing computation obtained as mentioned above.

In the present embodiment, an explanation is made in a case where an embedding degree is 12. However, also in a case where the embedding degree is 8 as described in claim 22 to claim 24, it is possible to compute Twisted Xate pairing at high speed by means of the similar algorithm.

An explanation is made in a case where, that is, the embedding degree is 8, and r(χ) and t (χ) are given as follows. r(χ)=9χ⁴+12χ³+8χ²+4χ+1 t(χ)=−9χ³−3χ²−2χ. In this case, representation of integer variable χ using p² and p³ with p as a characteristic is as below. p ³ =p ²+3χ+1

In this case, f′ _(χ,S)(Q) becomes, using l₅ which is a value at a rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (3χS, S), l₆ which is a value at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3χ+1)S), l₇ which is a value at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (3χp³S, p³S), and l₈ which is a value at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3χ+1)p³S), as follows.

[F78] f′ _(χ,S)(Q)=(f _(3χ,S)(Q)·l ₅ ·l ₆)^(p) ³ ·f _(3χ,p) ³ _(S)(Q)·l ₇ l ₈

In the similar way as in the case where the embedding degree is 12, in the case where the embedding degree is 8, an authentication server performs pairing computation by means of the started pairing computation program based on a flowchart shown in FIG. 4.

The electronic computer 10, causing the CPU 11 to function as an input unit by means of the started pairing computation program, inputs data of an integer variable χ and data of rational point S both of which are stored preliminarily in the memory unit 13, and data of rational point Q temporarily stored in the memory unit 13 as signature data into predetermined registers provided inside of the CPU 11 respectively. (Step T1).

Next, the electronic computer 10, causing the CPU 11 to function as a first computation unit by means of the pairing computation program, computes a rational function f_(3χ,S)(Q) by means of Miller's algorithm and stores a result of the computation in the memory unit 13. (Step T2). In the step T2, the electronic computer 10 performs the same computation processing as in a flowchart shown in FIG. 3, except that anti-logarithm of logarithm to base two instep S21 in the flowchart shown in FIG. 3 is 3χ and the first expression in step S22 is

[F79] λ_(T,T)←(3x_(T) ²+a)/2y_(T).   1. Here, “a” is a coefficient of first degree in an elliptic curve given as y²=x³+ax, aεF_(p).

Further, also in the step T2, the electronic computer 10 computes 3χS as well as the rational function f_(3χ,S)(Q) and stores results of the computation in the memory unit 13.

After computing rational function f_(3χ,S)(Q) and 3χS, the electronic computer 10, causing the CPU 11 to function as a second computation unit by means of the pairing computation program, reads a value of 3χS stored in the memory unit 13, performs computation of rational points p²(S) and (3χ+1)S respectively and stores values of the rational points in the memory unit 13 respectively. (Step T3).

Next, the electronic computer 10, causing the CPU 11 to function as a third computation unit by means of the pairing computation program, reads values of the respective rational points obtained in step T3 from the memory unit 13 respectively, computes l₅ which is a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χS, S) and l₆ which is a value at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3χ+1)S) respectively, and stores values of l₅ and l₆ in the memory unit 13 respectively. (Step T4).

In the same way as in the case where the embedding degree is 12, to be more specific, the electronic computer 10 computes l=1 _(A,B)(Q) as follows. λ_(A,B)←(y_(B)−y_(A))/(x_(B)−x_(A)), l_(A,B)(Q)←(x_(Q)−x_(B))λ_(A,B)−(y_(Q)−y_(B))

Next, the electronic computer 10, causing the CPU 11 to function as the second computation unit by means of the pairing computation program, computes rational point p³S (Step T5). To be more specific, the electronic computer reads values of rational points obtained in step T2 and step T3 from the memory unit 13, computes as below, and stores the value of p³S in the memory unit 13. p ³ S=p ² S +(3χ+1)S.

Next, the electronic computer 10, causing the CPU 11 to function as the first computation unit by means of the pairing computation program, computes a rational function f_(3χ,p^3S)(Q) by means of Miller's algorithm and stores a result of the computation in the memory unit 13. (Step T6).

The electronic computer 10 performs a computation of the rational function f_(3χ, p^3S)(Q) specifically as shown in FIG. 3. In particular, in step T6, the electronic computer 10 performs computation of 3χp³S as well as computation of the rational function f_(3χ,p^3S)(Q) and stores results of the computation in the memory unit 13.

Next, the electronic computer 10, causing the CPU 11 to function as the second computation unit by means of the pairing computation program, reads the value of 3χp³S stored in the memory unit 13, performs computation of rational points p²(p³S), (3χ+1)p³S respectively and stores values of the respective rational points in the memory unit 13 respectively. (Step T7).

Next, the electronic computer 10, causing the CPU 11 to function as the third computation unit by means of the pairing computation program, reads values of respective rational points obtained in step T7 from the memory unit 13 respectively, computes l₇ which is a value at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χp³S, p³S) and l₈ which is in a value at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3χ+1)p³S) respectively, and stores values of l₇ and l₈ in the memory unit 13 respectively. (Step T8).

Next, the electronic computer 10, causing the CPU 11 to function as a fourth computation unit by means of the pairing computation program, reads a result of the computation in the first computation unit (values of f_(3χ,S)(Q) and f_(3χ,p^3S)(Q)), and a result of the computation in the third computation unit (values of l₅, l₆, and l₇, l₈) from the memory unit 13 and computes f′ _(χ,S)(Q) expressed as below and stores a result of the computation in the memory unit 13 (Step T9).

[F80] f′ _(χ,S)(Q)=(f _(3χ,S)(Q)·l ₅ ·l ₆)^(p) ³ ·f _(3χp) ³ _(S)(Q)·l ₇ ·l ₈   (Equation 30)

To be more specific, the electronic computer 10 performs computation of f′ _(χ,S)(Q) as follows.

input : SεG₁, QεG₂, 3 χ, p output : f=f’ _(χ,) _(S)(Q) procedure :  1. A←f_(3χ,) _(S) //T₁←3 χ S is obtained  2. B←g_(T1, S)  3. T₂←T₁ +S  4. T₃←p²S  5. C←g_(T3, T2)  6. f←A·B·C  7. T₄←T₃ +T₂ //T₄←p³S  8. A←f_(3χ,) _(T4) ///T₁←3 χ p³S  9. B←g_(T1, T4) 10. T₂←T₁ +T4 11. T₃←p²T₄ 12. C←g_(T3, T2) 13. f←f

p³ ·A·B·C 14. Return f That is, the CPU 11 performs respective computations of the procedure described above using the memory unit 13, obtains a value of f′ _(χ,S)(Q) and stores the value in the memory unit 13. To be more specific, in step 1, the CPU 11 obtains f_(3χ,S) by means of Miller's algorithm, assigns the result to A, and also assigns 3χS obtained simultaneously to T₁. In step 2, the CPU 11 obtains l₅ which is a value at rational point Q of a straight line passing through rational points (3χS, S) and assigns the result to B. In step 3 and step 4, the CPU 11 obtains rational points (3χ+1)S and p²(S) and assigns the results to T₂ and T₃ respectively. In step 5, the CPU 11 obtains l₆ which is a value at rational point Q of a straight line passing through rational points (p²(S), (3χ+1)S) and assigns the result to C. In step 6, the CPU 11 takes a product of A, B and C, obtains a ( ) part of (Equation 30) described above and assigns the result to f. In step 7, the CPU 11 obtains p³S using relationship that p³=p²+3χ+1 and assigns the result to T₄. In step 8, the CPU 11 obtains f_(3χ,p^3S) by means of Miller's algorithm and assigns the result to A and also assigns 3χp³S obtained simultaneously to T₁. In step 9, the CPU 11 obtains l₇ which is a value at rational point Q of a straight line passing through rational points (3χp³S, p³S) and assigns the result to B. In step 10 and step 11, the CPU 11 obtains rational points (3χ+1)p³S, p²(p³S) and assigns the results to T₂ and T₃ respectively. In step 12, the CPU 11 obtains l₈ which is a value at rational point Q of a straight line passing through rational points (p²(p³S), (3χ+1)p³S) and assigns the result to C. In step 13, the CPU 11 obtains f′ _(χ,S)(Q) and assigns the result to f.

Next, the electronic computer 10, causing the CPU 11 to function as a fifth computation unit by means of the pairing computation program, reads a value of f′ _(χ,S)(Q) obtained in step T9 from the memory unit 13, performs an exponentiation of the final exponentiation in pairing e(S, Q) and stores a result of the computation in the memory unit 13. (Step T10).

Hereinafter, an embodiment of Twisted Xate multi-pairing in a case where an embedding degree is 12 is explained in detail. Explanation of the parts which are common to Twisted Xate pairing is omitted. Based on FIG. 5, Twisted Xate multi-pairing in a case where the embedding degree is 12 is explained.

The electronic computer 10, causing the CPU 11 to function as an input unit by means of the started pairing computation program, inputs data of an integer variable χ and data of rational point S both of which are stored preliminarily in the memory unit 13, and data of rational point Q temporarily stored in the memory unit 13 as signature data into predetermined registers provided inside of the CPU 11 respectively. (Step U1).

Next, the electronic computer 10, causing the CPU11 to function as a second computation unit by means of the pairing computation program, performs a computation of rational point pS and stores the value of pS in the memory unit 13 (step U2).

Next, the electronic computer 10, causing the CPU11 to function as the second computation unit by means of the pairing computation program, performs a computation of rational point pQ and stores the value of pQ in the memory unit 13 (step U3).

Next, the electronic computer 10, causing the CPU 11 to function as a first computation unit by means of the pairing computation program, computes a rational function F_(2χ,2S)(Z_(Q)) by means of Miller's algorithm with respect to multi-pairing and stores a result of the computation in the memory unit 13. (Step U4).

The electronic computer 10 performs computation of the rational function F_(2χ,2S) (Z_(Q)) specifically as shown in FIG. 6. In particular, in Step U4, the electronic computer 10 computes 2χS and 2χpS as well as the rational function F_(2χ,ZS)(Z_(Q)) and stores results of the computation in the memory unit 13.

That is, firstly, the electronic computer 10, as initial setting, sets f to be 1, T₁ to be S, T₂ to be pS and also sets the number of bits in a case of binary representation of the integer variable 2χ to be i (Step U41).

Next, the electronic computer 10 performs predetermined computation of the rational function F_(2χ,ZS) (Z_(Q)) part in step U42, performs predetermined computation of 2χS and 2χpS part in step U43, determines whether u_(i) which is a value of i-th bit counted from the least degree of the integer variable 2χ is “1” or “0” (Step U44) and, when u_(i)=1, further performs predetermined computation of F_(2χZS)(Z_(Q)) part in step U45 and predetermined computation of 2χS and 2χpS part in step U46.

Next, the electronic computer 10, when i≠1 (Step U47), decrements i (Step U48), returns to Step U42 and, by repeating the processing until i=1 in Step U47, computes the rational function F_(2χ,ZS)(Z_(Q), and 2χS and 2χpS, and stores results of the computation in the memory unit 13 respectively.

Next, the electronic computer 10, causing the CPU 11 to function as the second computation unit by means of the pairing computation program, reads values of 2χS and 2χpS stored in the memory unit 13, performs computation of rational points −S, (2χ−1)S, p¹⁰((2χ−1)S), 31 pS, (2χ−1)pS, and p¹⁰((2 χ−1)pS) respectively and stores values of the respective rational points in the memory unit 13 respectively. (Step U5).

Next, the electronic computer 10, causing the CPU 11 to function as a third computation unit by means of the pairing computation program, reads values of the respective rational points obtained in step U5 from the memory unit 13 respectively, computes l₁ which is a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)S, −S), l₂ which is a value at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2χ−1)pS, −pS), l₃ which is a value at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)S), 2χS) and l₄ which is a value at the rational point Q(x(_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2χ−1)pS), 2χpS) respectively, and stores values of l₁, l₂, l₃, and l₄ in the memory unit 13 respectively. (Step U6).

To be more specific, the electronic computer 10 computes l=l_(A,B)(Q) as follows. λ_(A,B)←(y_(B)−y_(A))/(x_(B)−x_(A)), l_(A,B)(Q)←(x_(Q)−x_(B))λ_(A,B)−(u_(Q)−y_(B))

That is, the CPU 11 computes a gradient λ_(A,B) using xy coordinates of points A and B, computes l_(A,B)(Q) using a result of the computation and xy coordinates of the points Q and B, and stores a result of the computation in the memory unit 13.

Next, the electronic computer 10, causing the CPU 11 to function as a fourth computation unit by means of the pairing computation program, reads a result of the computation in the first computation unit (Step U4), and a result of the computation in the third computation unit (values of l₁, l₂, l₃, and l₄) from the memory unit 13, computes f′ _(χ,S)(Q) expressed as below, and stores a result of the computation in the memory unit 13 (Step U7).

[F81] f′ _(χ,S)(Q)={F _(2χ,Z) _(S) (Z _(Q))·{l ₁ ^(p) ·l ₂}⁻¹}^(p) ¹⁰ ·F _(2χ,Z) _(S) (Z _(Q))·l ₃ ^(p) ·l ₄   (Equation 31)

Next, the electronic computer 10, causing the CPU 11 to function as a fifth computation unit by means of the pairing computation program, reads a value of f′ _(χ,S)(Q) obtained in step U7 from the memory unit 13, performs an exponentiation of the final exponentiation in pairing e(S, Q), and stores a aresult of the computation in the memory unit 13 (Step U8).

To be more specific, the electronic computer 10 performs a computation of f′ _(χ,S)(Q) as follows.

input : Z_(s)={S, pSεG₁}, Z_(Q)={pQ, QεG₂}, 2 χ, p output : f=f’ _(χ,) _(S)(Q) procedure :  1. A←MMA (2 χ, 2, Z_(S), Z_(Q) ) //T₁←2 χ S, T₂←2 χ pS  2. T₃← −S  3. T₄←T₁ +T₃  4. B←g_(T4, T3)  5. T₅← −pS  6. T₆←T₂ +T₅  7. C←g_(T5, T6)  8. f←B^(p) ·C  9. f←A ·f⁻¹ 10. T₃←p¹⁰T₄ 11. B←g_(T3, T1) 12. T₅←p¹⁰T₆ 13. C←g_(T5, T2) 14. D←B^(p) ·C 15. f←f

p¹⁰ ·A ·D 16. Return f That is, the CPU 11 performs respective computations of the procedure described above using the memory unit 13, obtains f′ _(χ,S)(Q), and stores the value in the memory unit 13.

To be more specific, in step 1, the CPU 11 obtains F_(2χ,ZS) by means of Miller's algorithm with respect to multi-pairing, assigns the result to A, and also assigns 2χS and 2χpS obtained simultaneously to T₁ and T₂ respectively. In step 2 and step 3, the CPU 11 obtains rational points −S and (2χ1)S and assigns the result to T₃ and T₄ respectively. In step 4, the CPU 11 obtains l₁ which is a value at rational point Q of a straight line passing through rational points ((2χ−1)S, −S) and assigns the result to B. In step 5 and 6, the CPU 11 obtains rational points −pS and (2χ−1)pS and assign the result to T₅ and T₆ respectively. In step 7, the CPU 11 obtains l₂ which is a value at rational point Q of a straight line passing through rational points (−pS, (2χ−1)pS) and assigns the result to C. In step 8, the CPU 11 takes a product of B to the power of p and C, obtains a {} part of the inside of (Equation 31) described above and assigns the result to f. In step 9, the CPU 11 obtains a {} part of the outside of (Equation 31) described above and assigns the result to f. In step 10, the CPU 11 obtains a value of rational point p¹⁰((2 χ−1)S) and assigns the result to T₃. In step 11, the CPU 11 obtains l₃ which is a value at rational point Q of a straight line passing through rational points (p¹⁰((2χ−1)S), 2χS) and assigns the results to B. In step 12, the CPU 11 obtains rational point p¹⁰((2χ−1)pS) and assigns the result to T₅. In step 13, the CPU 11 obtains l₄ which is a value at rational point Q of a straight line passing through rational points (p¹⁰((2χ−1)pS), 2χpS) and assigns the result to C. In step 14, the CPU 11 obtains a product of B to the power of p and C. In step 15, the CPU 11 obtains the whole of (Equation 31) described above and assigns the result (f′ _(χ,S)(Q)) to f.

In the embodiment, Twisted Xate multi-pairing in the case where the embedding degree is 12 is explained. However, it is possible to compute Twisted Xate multi-pairing at high speed also in a case where the embedding degree is 8 by means of the similar algorithm. Explanation of the parts which are common to Twisted Xate pairing is omitted. Computation of multi-pairing in the case where the embedding degree is 8 is explained with reference to FIG. 7.

The electronic computer 10, causing the CPU 11 to function as an input unit by means of the started pairing computation program, inputs data of an integer variable χ and data of rational point S both of which are stored preliminarily in the memory unit 13, and data of rational point Q temporarily stored in the memory unit 13 as signature data into predetermined registers provided inside of the CPU 11 respectively. (Step V1).

Next, the electronic computer 10, causing the CPU 11 to function as a second computation unit by means of the pairing computation program, performs a computation of rational point p³S and stores the value of p³S in the memory unit 13 (step V2).

Next, the electronic computer 10, causing the CPU 11 to function as the second computation unit by means of the pairing computation program, performs a computation of rational point p³Q and stores the value of p³Q in the memory unit 13 (step V3).

Next, the electronic computer 10, causing the CPU 11 to function as a first computation unit by means of the pairing computation program, computes a rational function F_(3χ,ZS)(Z_(Q)) by means of Miller's algorithm with respect to multi-pairing and stores a result of the computation in the memory unit 13. (Step V4). In this Step V4, the electronic computer 10 performs the same computation processing as a flowchart shown in FIG. 6 except that anti-logarithm of logarithm to base two is 3χ, and all p become p³.

Also in step V4, the electronic computer 10 computes 3χS and 3χp³S as well as the rational function F_(3χ,ZS)(Z_(Q)) and stores results of the computation in the memory unit 13.

After computing the rational function F_(3χ,ZS) (Z_(Q)), and 3χS and 3χp³S, the electronic computer 10, causing the CPU 11 to function as the second computation unit by means of the pairing computation program, reads values of 3χS and 3χp³S stored in the memory unit 13, performs computation of rational points p²(S), (3χ+1)S, p²(p³S), and (3χ+1)p³S respectively and stores values of the respective rational points in the memory unit 13 respectively. (Step V5).

Next, the electronic computer 10, causing the CPU 11 to function as a third computation unit by means of the pairing computation program, reads values of the respective rational points obtained in step V5 from the memory unit 13 respectively, computes l₅ which is a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3χS, S), l₆ which is a value at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3χ+1)S), l₇ which is a value at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (3χp³S, p³S) and l₈ which is a value at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3χ+1)p³S) respectively, and stores values of l₅, l₆, l₇, and l₈ in the memory unit 13 respectively. (Step V6).

To be more specific, in the same way as in the case where the embedding degree is 12, the electronic computer 10 computes l=l_(A,B)(Q) as follows. λ^(A,B)←(y_(B)−y_(A))/(x_(B)-X_(A)), l_(A,B)(Q)←(x_(Q)−x_(B))λ_(A,B)−(y_(Q)−y_(B))

Next, the electronic computer 10, causing the CPU 11 to function as a fourth computation unit by means of the pairing computation program, reads a result of the computation in the first computation unit, and a result of the computation in the third computation unit (values of l₅, l₆, l₇, and l₈) from the memory unit 13 and computes f′ _(χ,S)(Q) expressed as below and stores a result of the computation in the memory unit 13 (Step V7).

[F82] f′ _(χ,S)(Q)=F _(3χ,Z) _(S) ·{l ₅ ·l ₆}^(p) ^(s) ·l ₇ ·l ₈   (Equation 32)

To be more specific, the electric computer 10 performs the computation of f′ _(χ,S)(Q) as follows.

input : Z_(S)={S, p³SεG₁}, Z_(Q)={p³Q, QεG₂}, 3 χ, p output : f=f’ _(χ,) _(S)(Q) procedure :  1. A←MMA (3 χ, 2, Z_(S), Z_(Q) ) //T₁←3 χ S, T₂←3 χ p³S  2. B←g_(3 χS, S)  3. T₃←T₁ +S  4. T₄←p²S  5. C← g_(T4, T3)  6. f←B ·C  7. f←f 

p³  8. f←A ·f  9. B← 

10. T₃←T₂ +p³S 11. T4←p²p³S 12. C←g_(T4, T3) 13. f←f · B ·C 14. Return f That is, the CPU 11 performs respective computations of the procedure described above using the memory unit 13, obtains f′ _(χ,S)(Q), and stores the value in the memory unit 13. To be more specific, In step 1, the CPU 11 obtains F_(3χ,ZS) by means of Miller's algorithm with respect to multi-pairing and assigns the result to A and also assigns 3χS and 3χp³S obtained simultaneously to T₁ and T₂ respectively. In step 2, the CPU 11 obtains l₅ which is a value at rational point Q of a straight line passing through rational points (3χS, S) and assigns the result to B. In step 3 and step 4, the CPU 11 obtains rational points (3 χ+1)S and p²(S) and assigns the result to T₃ and T₄ respectively. In step 5, the CPU 11 obtains l₆ which is a value at rational point Q of a straight line passing through rational points (p²(S), (3χ+1)S) and assigns the result to C. In step 6, the CPU 11 takes a product of B and C, obtains an inside of { } part of (Equation 32) described above and assigns the result to f. In step 7, the CPU 11 obtains a {} part of (Equation 32) described above. In step 8, the CPU 11 obtains a result up to the {} part of (Equation 32) described above and assign the result to f. In step 9, the CPU 11 obtains l₇ which is a value at rational point Q of a straight line passing through rational points (3χp³S, p³S) and assigns the result to B. In step 10 and step 11, the CPU 11 obtains rational points (3χ+1) p³S and p²(p³S), and assigns the results to T₃ and T₄ respectively. In step 12, the CPU 11 obtains l₈ which is a value at rational point Q of a straight line passing through rational points (p²(p³S), (3χ+1)p³S) and assigns the result to C. In step 13, the CPU 11 obtains f′ _(χ,S)(Q), and assigns the result to f.

Next, the electronic computer 10, causing the CPU 11 to function as a fifth computation unit by means of the pairing computation program, reads the value of f′ _(χ,S)(Q) obtained in step V7, performs an exponentiation of the final exponentiation in pairing e(S, Q), and stores a result of the computation in the memory unit 13 (Step V8).

In addition, in the explanation made above, the memory unit 13 is used as storage means for storing each computation result. However, registers provided with inside of the CPU 11 may be used as the storage means described above.

EXPLANATION OF SYMBOLS

10 electronic computer

11 CPU

12 storage unit

13 memory unit

14 bus

15 input output control part

20 telecommunication line

30 client device 

What is claimed is:
 1. A pairing computation device, wherein an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), the pairing computation device comprising a CPU which computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and F be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA), wherein the order r and a trace t of the Frobenius endomorphism φ_(p) are specified preliminarily according to the embedding degree k using the integer variable χ, and the CPU includes: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a computation unit which computes F; a computation unit which computes a value of a straight line passing through given rational points at a rational point Q (x_(Q), y_(Q)); a computation unit which computes f′ _(χ,) ^(S) (Q) using said F and said value ; and a computation unit which computes the pairing e(S,Q) using said f′ _(χ,) ^(S) (Q) as [F1] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ^(k) ^(−1)/r).
 2. A pairing computation method, wherein an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), an electronic computer which includes a CPU computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and F be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA), the order r and a trace t of the Frobenius endomorphism φ_(p) being specified preliminarily according to the embedding degree k using the integer variable χ, the pairing computation method comprising: a step of inputting the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit; a step of computing F by causing the CPU of the electronic computer to function as a computation unit; a step of computing a value of a straight line passing through given rational points at a rational point Q(x_(Q), y_(Q)) by causing the CPU of the electronic computer to function as a computation unit; a step of computing f′ _(χ,) ^(S) (Q) using said F and said value by causing the CPU of the electronic computer to function as a computation unit; and a step of computing the pairing e(S,Q) using said f′ _(χ,) ^(S) (Q) as [F2] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ^(k) _(−1)/r) by causing the CPU of the electronic computer to function as a computation unit.
 3. A non-transitory computer recording medium storing a pairing computation program, wherein an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and F be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA), the order r and a trace t of the Frobenius endomorphism φ_(p) being specified preliminarily according to the embedding degree k using the integer variable χ, the pairing computation program causing the CPU of the electronic computer to function as: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a computation unit which computes F; a computation unit which computes a value of a straight line passing through given rational points at a rational point Q(x_(Q), y_(Q)); a computation unit which computes f′ _(χ,) ^(S) (Q) using said F and said value; and a computation unit which computes the pairing e(S,Q) using said f′ _(χ,) ^(S) (Q) as [F3] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ^(k) ^(−1)/r).
 4. A pairing computation device, wherein an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ¹²/(F* _(p) ¹²)^(r), the pairing computation device comprising a CPU which computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and pS, Z_(Q) be a set of rational points pQ and Q, and F_(2χ,ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA), wherein the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=36_(χ) ⁴−36_(χ) ³+18_(χ) ²−6_(χ)+1, t(χ)=6_(χ) ²+1, and a representation of the integer variable χ using p¹⁰ with p as a characteristic is p≡(2_(χ)−1)p ¹⁰+2χ(modr(χ)) and the CPU includes: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a first computation unit which computes F_(2χ,ZS)(Z_(Q)); a second computation unit which computes given rational points using 2_(χ)S and 2_(χ)pS which are calculated when computing said F_(2χ,ZS)(Z_(Q)); a third computation unit which computes a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through the given rational points; a fourth computation unit which computes f′ _(χ,) ^(S) (Q) using said F_(2χ,ZS)(Z_(Q)) and said value; and a fifth computation unit which computes the pairing e(S,Q) using said f′ _(χ,) ^(S) (Q) as [F4] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ¹² ^(−1)/r).
 5. The pairing computation device according to claim 4, wherein the second computation unit computes respective rational points −S, (2_(χ)−1)S, p¹⁰((2_(χ)−1)S), −pS, (2_(χ)−1)pS, p¹⁰((2_(χ)−1)pS) in order using previously obtained results, the third computation unit respectively computes a value l₁ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)S, −S), a value l₂ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)pS, −pS), a value l₃ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1)S), 2_(χ)S), and a value l₄ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1)pS), 2_(χ)pS), and the fourth computation unit computes f′ _(χ,) ^(S) (Q) as [F5] f′ _(χ,) ^(S) (Q)={F ₂ _(χ) _(,Z) _(S) (Z _(Q))·{l ₁ ^(p) ·l ₂}⁻¹}^(p) ¹⁰ ·F ₂ _(χ) _(,Z) _(S) (Z _(Q))·l ₃ ^(p) ·l ₄.
 6. A pairing computation method, wherein an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂ →F* _(p) ¹²/(F* _(p) ¹²)^(r), an electronic computer which includes a CPU computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and pS, Z_(Q) be a set of rational points pQ and Q, and F_(2χ,ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA), the order r and a trace t of the Frobenius endomorphism φ_(p) being specified using the integer variable χ as, r(χ)=36_(χ) ⁴−36_(χ) ³+18_(χ) ²−6_(χ)+1 t(χ)=6_(χ) ²+1, and a representation of the integer variable χ using p¹⁰ with p as a characteristic being p≡(2_(χ)−1)p ¹⁰+2_(χ)(modr(χ)), the pairing computation method comprising: an input step which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit; a first computation step which computes F_(2χ,ZS)(Z_(Q)) by causing the CPU of the electronic computer to function as a first computation unit; a second computation step which computes given rational points using 2χS and 2χpS which are calculated when computing said F_(2χ,ZS)(Z_(Q)) by causing the CPU of the electronic computer to function as a second computation unit; a third computation step which computes a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through the given rational points by causing the CPU of the electronic computer to function as a third computation unit; a fourth computation step which computes f′_(χ,) ^(S) (Q) using said F_(2χ, ZS)(Z_(Q)) and said value by causing the CPU of the electronic computer to function as a fourth computation unit ; and a fifth computation step which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F6] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ¹² ^(−1)/r), by causing the CPU of the electronic computer to function as a fifth computation unit.
 7. The pairing computation method according to claim 6, wherein the second computation step computes respective rational points −S, (2_(χ)−1)S, p¹⁰((2_(χ)−1)S), −pS, (2_(χ)−1)pS, p¹⁰((2_(χ)−1)pS) in order using previously obtained results, the third computation step respectively computes a value l₁ at a rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)S, −S), a value l₂ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)pS, −pS), a value l₃ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰ ((2_(χ)−1) S), 2_(χ)S), and a value l₄ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1)pS), 2_(χ)pS), and the fourth computation step computes f′_(χ,) ^(S) (Q) as [F7] f′ _(χ,) ^(S) (Q)={F ₂ _(χ) _(,Z) _(S) (Z _(Q))·{l ₁ ^(p) ·l ₂}⁻¹}^(p) ¹⁰ ·F ₂ _(χ) _(,Z) _(S) (Z _(Q))·l ₃ ^(p) ·l ₄.
 8. A non-transitory computer recording medium storing a pairing computation program, wherein an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ¹²/(F* _(p) ¹²)^(r), the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and pS, Z_(Q) be a set of rational points pQ and Q, and F_(2χ,ZS) (Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA), the order r and a trace t of the Frobenius endomorphism φ_(p) being specified using the integer variable χ as, r(χ)=36_(χ) ⁴−36_(χ) ³+18_(χ) ²−6_(χ)+1, t(χ)=6_(χ) ²+1, and a representation of the integer variable χ using p¹⁰ with p as a characteristic being p≡(2_(χ)−1)p¹⁰+2_(χ)(modr(χ)), the pairing computation program causing the CPU of the electronic computer to function as: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a first computation unit which computes F_(2χ,ZS)(Z_(Q)); a second computation unit which computes given rational points using 2_(χ)S and 2_(χ)pS which are calculated when computing said F_(2χ,ZS)(Z_(Q)); a third computation unit which computes a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through the given rational points; a fourth computation unit which computes f′_(χ,) ^(S) (Q) using said F_(2χ,ZS)(Z_(Q)) and said value; and a fifth computation unit which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F8] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ¹² ^(−1)/r).
 9. The non-transitory computer recording medium storing a pairing computation program according to claim 8, wherein the pairing computation program causes: the CPU of the electronic computer which functions as the second computation unit to compute respective rational points −S, (2_(χ)−1)S, p¹⁰((2_(χ)−1)S), −pS, (2_(χ)−1)pS, p¹⁰((2_(χ)−1)pS) in order using previously obtained results; the CPU of the electronic computer which functions as the third computation unit to compute respectively a value l₁ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)S, −S), a value l₂ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)pS, −pS), a value l₃ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1)S), 2_(χ)S), and a value l₄ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1)pS), 2_(χ)pS); and the CPU of the electronic computer which functions as the fourth computation unit to compute f′_(χ,) ^(S) (Q) as [F9] f′ _(χ,) ^(S) (Q)={F ₂ _(χ,) _(Z) _(S) (Z _(Q))·{l ₁ ^(p) ·l ₂}⁻¹}^(p) ¹⁰ ·F ₂ _(χ,) _(Z) _(S) (Z _(Q))·l ₃ ^(p) ·l ₄.
 10. A pairing computation device, wherein an elliptic curve is given as y²=x³+ax, aεF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ⁸/(F* _(p) ⁸)^(r), the pairing computation device comprising a CPU which computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and p³S, Z_(Q) be a set of rational points p³Q and Q, and F_(3χ,ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA), wherein the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=9_(χ) ⁴+12_(χ) ³+8_(χ) ²+4_(χ)+1, t(χ)=9_(χ) ³−3_(χ) ²−2_(χ), and a representation of the integer variable χ using p² and p³ with p as a characteristic is p ³≡p²+3_(χ)+1(modr(χ)) and the CPU includes: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a first computation unit which computes F_(3 χ,ZS)(Z_(Q)); a second computation unit which computes respective rational points p²(S), p²(p³S), (3_(χ)+1)S, (3_(χ)+1)p³S in order using previously obtained results; a third computation unit which respectively computes a value l₅ at a rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (3_(χ)S, S), a value l₆ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3_(χ)+1)S), a value l₇ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (3_(χp) ³S, p³S), and a value l₈ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3_(χ)+1)p³S); and the fourth computation unit which computes f′_(χ,) ^(S) (Q) using a computation result of the first computation unit and a computation result of the third computation unit as [F10] f′ _(χ,) ^(S) (Q)=F ₃ _(χ,) _(Z) _(S) (Z _(Q)){l ₅ ·l ₆}^(p) ³ ·l ₇ ·l ₈; and a fifth computation unit which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F11] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ³ ^(−1)/r).
 11. A pairing computation method, wherein an elliptic curve is given as y²=x³+ax, aεF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ⁸/(F* _(p) ⁸)^(r), an electronic computer which includes a CPU computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and p³S, Z_(Q) be a set of rational points p³Q and Q, and F_(3χ,ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA), the order r and a trace t of the Frobenius endomorphism φ_(p) being specified using the integer variable χ as, r(χ)=9_(χ) ⁴+12_(χ) ³+8_(χ) ²+4_(χ)+1. t(χ)=−9_(χ) ³−3_(χ) ²−2_(χ), and a representation of the integer variable χ using p² and p³ with p as a characteristic being p ³≡p²+3_(χ)+1(modr(χ)) the pairing computation method comprising: an input step which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit; a first computation step which computes F_(3χ,ZS)(Z_(Q)) by causing the CPU of the electronic computer to function as a first computation unit; a second computation step which computes respective rational points p²(S), p²(p³S), (3_(χ)+1)S, (3_(χ)+1)p³S in order using previously obtained results by causing the CPU of the electronic computer to function as a second computation unit; a third computation step which respectively computes a value l₅ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3_(χ)S, S), a value l₆ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3_(χ)+1)S), a value l₇ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3_(χp) ³S, p³S), and a value l₈ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3_(χ)+1)p³S) by causing the CPU of the electronic computer to function as a third computation unit; and a fourth computation step which computes f′_(χ,) ^(S) (Q) using said F_(3χ,ZS)(Z_(Q)) and said values l₅, l₆, l₇, l₈ as [F12] f′ _(χ,) ^(S) (Q)=F ₃ _(χ,) _(Z) _(S) (Z _(Q)){l ₅ ·l ₆}^(p) ³ ·l ₇ ·l ₈ by causing the CPU of the electronic computer to function as a fourth computation unit; and a fifth computation step which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F13] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ³ ^(−1)/r) by causing the CPU of the electronic computer to function as a fifth computation unit.
 12. A non-transitory computer recording medium storing a pairing compututation program, wherein an elliptic curve is given as y²=x²+ax, aεF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ⁸/(F* _(p) ⁸)^(r), the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, Z_(S) be a set of rational points S and p³S, Z_(Q) be a set of rational points p³Q and Q, and F_(3χ,ZS)(Z_(Q)) be a rational function which is calculated using Miller's algorithm with respect to multi-pairing (MMA), the order r and a trace t of the Frobenius endomorphism φ_(p) being specified using the integer variable χ as, r(χ)=9_(χ) ⁴+12_(χ) ³+8_(χ) ²+4_(χ)+1, t(χ)=−9_(χ) ³−3_(χ) ²−2_(χ), and a representation of the integer variable χ using p² and p³ with p as a characteristic being p ³≡p²+3_(χ+1)(modr(χ)) the pairing computation program causing the CPU of the electronic computer to function as: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a first computation unit which computes F_(3χ,ZS)(Z_(Q)); a second computation unit which computes respective rational points p²(S), p²(p³S), (3_(χ)+1)S, (3_(χ)+1)p²S in order using previously obtained results; a third computation unit which respectively computes a value l₅ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3_(χ)S, S), a value l₆ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3χ+1)S), a value l₇ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3_(χp) ³S, p³S), and a value l₈ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3_(χ)+1)p³S); and a fourth computation unit which computes f′_(χ,) ^(S) (Q) using said F_(3χ,ZS)(Z_(Q)) and said values l₅, l₆, l₇, l₈ as [F14] f′ _(χ,) ^(S) (Q)=F ₃ _(χ,) _(Z) _(S) (Z _(Q)){l ₅ ·l ₆}^(p) ³ ·l ₇ · ₈; and a fifth computation unit which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F15] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ³ ^(−1)/r).
 13. A pairing computation device, wherein an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), the pairing computation device comprising a CPU which computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f be a rational function which is calculated using Miller's algorithm, wherein the order r and a trace t of the Frobenius endomorphism φ_(p) are specified preliminarily according to the embedding degree k using the integer variable χ, and the CPU includes: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a computation unit which computes f; a computation unit which computes a value of a straight line passing through given rational points at a rational point Q(x_(Q), y_(Q)); a computation unit which computes f′_(χ,) ^(S) (Q) using said f and said value; and a computation unit which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F16] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ^(k) ^(−1)/r).
 14. A pairing computation method, wherein an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), an electronic computer which includes a CPU computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f be a rational function which is calculated using Miller's algorithm, the order r and a trace t of the Frobenius endomorphism φ_(p) being specified preliminarily according to the embedding degree k using the integer variable χ, the pairing computation method comprising: a step of inputting the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit; a step of computing f by causing the CPU of the electronic computer to function as a computation unit; a step of computing a value of a straight line passing through given rational points at a rational point Q(x_(Q), y_(Q)) by causing the CPU of the electronic computer to function as a computation unit; a step of computing f′_(χ,) ^(S) (Q) using said f and said value by causing the CPU of the electronic computer to function as a computation unit; and a step of computing the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F17] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ^(k) ^(−1)/r) by causing the CPU of the electronic computer to function as a computation unit.
 15. A non-transitory computer recording medium storing a pairing computation program, wherein an elliptic curve is given as y²=x³+ax+b, aεF_(p), bεF_(p), letting k be an embedding degree, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ^(k), E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as non-degenerate bilinear map e:G ₁ ×G ₂ →F* _(p) ^(k)/(F* _(p) ^(k))^(r), the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f be a rational function which is calculated using Miller's algorithm, the order r and a trace t of the Frobenius endomorphism φ_(p) being specified preliminarily according to the embedding degree k using the integer variable χ, the pairing computation program causing the CPU of the electronic computer to function as: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a computation unit which computes f; a computation unit which computes a value of a straight line passing through given rational points at a rational point Q(x_(Q), y_(Q); a computation unit which computes f′_(χ,) ^(S) (Q) using said f and said value; and a computation unit which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F18] e(S,Q)=f′ _(χ,) ^(S) (Q)^(p) ^(k) ^(−1)/r).
 16. A pairing computation device, wherein an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ¹²/(F* _(p) ¹²)^(r), the pairing computation device comprising a CPU which computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, f₂ _(χ,) _(S)(Q) be a rational function which is calculated using Miller's algorithm, wherein the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=36_(χ) ⁴−36_(χ) ³+18_(χ) ²−6_(χ)+1, t(χ)=6_(χ) ²+1, and a representation of the integer variable χ using p¹⁰ with p as a characteristic is p≡(2χ−1)p ¹⁰+2χ(modr(χ)) and the CPU includes: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a first computation unit which computes f_(2χ,S)(Q) and f_(2χ,pS)(Q); a second computation unit which computes given rational points using 2_(χ)S and 2_(χ)pS which are calculated when computing said f_(2χ,S)(Q) and f_(2χ,pS)(Q); a third computation unit which computes a value at a rational point Q (x_(Q), y_(Q)) of a straight line passing through the given rational points; a fourth computation unit which computes f′_(χ,) ^(S) (Q) using said f_(2χ,S)(Q) and said value; and a fifth computation unit which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F19] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ¹² ^(−1)/r).
 17. The pairing computation device according to claim 16, wherein the second computation unit computes respective rational points −S, (2_(χ)−)S, p¹⁰((2_(χ)−1)S), −pS, (2_(χ)−1)pS, p¹⁰((2_(χ)−1)pS) in order using previously obtained results, the third computation unit respectively computes a value l₁ at a rational point Q (x_(Q), y_(Q)) on a straight line passing through rational points ((2_(χ)−1)S, −S), a value l₂ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1)S), 2_(χ)S), a value l₃ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)pS, −pS), and a value l₄ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1) pS), 2_(χ)pS), and the fourth computation unit computes f′_(χ,) ^(S) (Q) using the values of the rational point Q (x_(Q), y_(Q)) l₁, l₂, l₃, l₄ as [F20] f′ _(χ,) ^(S) (Q)=({f _(2χ) ^(,S) (Q)·l ₁ ⁻¹}^(p) ¹⁰ ·f _(2χ) ^(,S) (Q)·l ₂)^(p) ·{f _(2χ,pS)(Q)·l ₃ ⁻¹}^(p) ¹⁰ ·f _(2χ,pS)(Q)·l ₄.
 18. A pairing computation method, wherein an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ¹²/(F* _(p) ¹²)^(r), an electronic computer which includes a CPU computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f_(2χ) ^(,S) (Q) be a rational function which is calculated using Miller's algorithm, the order r and a trace t of the Frobenius endomorphism φ_(p) being specified using the integer variable χ as, r(χ)=36_(χ) ⁴−36_(χ) ³+18_(χ) ²−6_(χ)+1, t(χ)=6_(χ) ²+1, and a representation of the integer variable χ using p¹⁰ with p as a characteristic being p≡(2_(χ)−1)p ¹⁰+2_(χ)(modr(χ)), the pairing computation method comprising: an input step which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit; a first computation step which computes f_(2χ,) ^(S) (Q) and f_(2χ,pS)(Q) by causing the CPU of the electronic computer to function as a first computation unit; a second computation step which computes given rational points using 2_(χ)S and 2_(χ)pS which are calculated when computing said f_(2χ,) ^(S) (Q) and f_(2χ,pS)(Q) by causing the CPU of the electronic computer to function as a second computation unit; a third computation step which computes a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through the given rational points by causing the CPU of the electronic computer to function as a third computation unit; a fourth computation step which computes f′_(χ,) ^(S) (Q) using said f_(2χ,) ^(S) (Q) and said value by causing the CPU of the electronic computer to function as a fourth computation unit; and a fifth computation step which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F21] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ¹² ^(−1)/r). by causing the CPU of the electronic computer to function as a fifth computation unit.
 19. The pairing computation method according to claim 18, wherein the second computation step computes respective rational points −S, (2_(χ)−1)S, p¹⁰((2_(χ)−1)S), −pS, (2_(χ)−1)pS, p¹⁰((2_(χ)−1)pS) in order using previously obtained results, the third computation step respectively computes a value l₁ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)S, −S), a value l₂ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1)S), 2_(χ)S), a value l₃ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)pS, −pS), and a value l₄ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1)pS), 2_(χ)pS), and the fourth computation step computes f′_(χ,) ^(S) (Q) using the values of the rational point Q(x_(Q), y_(Q)) l₁, l₂, l₃, l₄ as [F22] f′ _(χ,) ^(S) (Q)=({f _(2χ,) ^(S) (Q)·l ₁ ⁻¹}^(p) ¹⁰ ·f _(2χ,) ^(S) (Q)·l ₂)^(p) ·{f _(2χ,pS)(Q)·l ₃ ⁻¹}^(p) ¹⁰ ·f _(2χ,pS)(Q)·l ₄.
 20. A non-transitory computer recording medium storing a pairing computation program, wherein an elliptic curve is given as y²=x³+b, bεF_(p), letting an embedding degree be 12, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ¹², E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ¹²/(F* _(p) ¹²)^(r), the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f_(2χ,) ^(S) (Q) and f_(2χ,pS)(Q) be a rational function which is calculated using Miller's algorithm, the order r and a trace t of the Frobenius endomorphism φ_(p) being specified using the integer variable χ as, r(χ)=36_(χ) ⁴−36_(χ) ³+18_(χ) ²−6_(χ)+1, t(χ)=6_(χ) ²+1, and a representation of the integer variable χ using p¹⁰ with p as a characteristic being p≡(2_(χ−1))p ¹⁰+2_(χ)(modr(χ)), the pairing computation program causing the CPU of the electronic computer to function as: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a first computation unit which computes f_(2χ,) ^(S) (Q) and f_(2χ,pS)(Q); a second computation unit which computes given rational points using 2_(χ)S and 2_(χ)pS which are calculated when computing said f_(2χ,) ^(S) (Q) and f_(2χ,pS)(Q); a third computation unit which computes a value at a rational point Q(x_(Q), y_(Q)) of a straight line passing through the given rational points; a fourth computation unit which computes f′_(χ,) ^(S) (Q) using said f_(2χ,) ^(S) (Q) and said value; and a fifth computation unit which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F23] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ¹² ^(−1)/r).
 21. The non-transitory computer recording medium storing a pairing computation program according to claim 20, wherein The pairing computation program causes: the CPU of the electronic computer which functions as the second computation unit to compute respective rational points −S, (2_(χ−)1)S, p¹⁰((2_(χ)−1)S), −pS, (2_(χ)−1)pS, p¹⁰((2_(χ)−1)pS) in order using previously obtained results; the CPU of the electronic computer which functions as the third computation unit to respectively compute a value l₁ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)S, −S), a value l₂ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1)S), 2_(χ)S), a value l₃ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points ((2_(χ)−1)pS, −pS), and a value l₄ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p¹⁰((2_(χ)−1)pS), 2_(χ)pS); and the CPU of the electronic computer which functions as the fourth computation unit to compute f′_(χ,) ^(S) (Q) using the values of the rational point Q (x_(Q), y_(Q)) l₁, l₂, l₃, l₄ as [F24] f′ _(χ,) ^(S) (Q)=({f _(2χ,) ^(S) (Q)·l ₁ ⁻¹}^(p) ¹⁰ ·f _(2χ,) ^(S) (Q)·l ₂)^(p) ·{f _(2χ,pS)(Q)·l ₃ ⁻¹}^(p) ¹⁰ ·f _(2χ,pS)(Q)·l ₄.
 22. A pairing computation device, wherein an elliptic curve is given as y²=x³+ax, aεF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]χKer(φ_(p)−[1]), G ₂ =E[r]χKer(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ⁸/(F* _(p) ⁸)^(r), the pairing computation device comprising a CPU which computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f_(3χ,) ^(S) (Q) be a rational function which is calculated using Miller's algorithm, wherein the order r and a trace t of the Frobenius endomorphism φ_(p) are specified using the integer variable χ as, r(χ)=9_(χ) ⁴+12_(χ) ³+8_(χ) ²+4_(χ)+1, t(χ)=−9_(χ) ³−3_(χ) ²−2_(χ), and a representation of the integer variable χ using p² and p³ with p as a characteristic is p ³≡p²+3_(χ)+1(modr(χ)) and the CPU includes: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a first computation unit which computes f_(3χ,) ^(S) (Q) and f_(3χ,p^3S)(Q); a second computation unit which computes respective rational points p²(S), (3_(χ)+1)S, p²(p³S), (3_(χ)+1)p³S in order using previously obtained results; a third computation unit respectively computes a value l₅ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3_(χ)S, S), a value l₆ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3_(χ)+1)S), a value l₇ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3_(χ)p³S, p³S), and a value l₈ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3_(χ)+1)p³S); and the fourth computation unit which computes f′_(χ,) ^(S) (Q) using the computation result of the first computation unit and the computation result of the third computation unit as [F25] f′ _(χ,) ^(S) (Q)=(f _(3χ,) ^(S) (Q)·l ₅ ·l ₆)^(p) ³ f _(3χ,p) ³ _(S)(Q)·l ₇ ·l ₈; and a fifth computation unit which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F26] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ³ ^(−1)/r).
 23. A pairing computation method, wherein an elliptic curve is given as y²=x³+ax, aεF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ⁸/(F* _(p) ⁸)^(r), an electronic computer which includes a CPU computes the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f_(3χ,) ^(S) (Q) be a rational function which is calculated using Miller's algorithm, the order r and a trace t of the Frobenius endomorphism φ_(p) being specified using the integer variable χ as, r(χ)=9_(χ) ⁴+12_(χ) ³+8_(χ) ²+4_(χ+1,) t(χ)=−9_(χ) ³−3_(χ) ²−2_(χ), and a representation of the integer variable χ using p² and p³ with p as a characteristic being p ³≡p²+3_(χ)+1(modr(χ)) the pairing computation method comprising: an input step which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers by causing the CPU of the electronic computer to function as an input unit; a first computation step which computes f_(3χ,) ^(S) (Q) and f_(3χ,p^3S)(Q) by causing the CPU of the electronic computer to function as a first computation unit; a second computation step which computes respective rational points p²(S), (3_(χ)+1)S, p²(p³S), (3_(χ)+1) p³S in order using previously obtained results by causing the CPU of the electronic computer to function as a second computation unit; a third computation step which respectively computes a value l₅ at a rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3_(χ)S, S), a value l₆ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3_(χ)+1)S), a value l₇ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (3_(χ)p³S, p³S), and a value l₈ at the rational point Q(x_(Q), y_(Q)) of a straight line passing through rational points (p²(p³S), (3_(χ)+1)p³S) by causing the CPU of the electronic computer to function as a third computation unit; and a fourth computation step which computes f′_(χ,) ^(S) (Q) using said f_(3χ,) ^(S) (Q) and said f_(3χ,p^3S)(Q) and said values l₅, l₆, l₇, l₈ as [F27] f′ _(χ,) ^(S) (Q)=(f _(3χ,) ^(S) (Q)·l ₅ ·l ₆)^(p) ³ f _(3χ,p) ₃ _(S)(Q)·l ₇ ·l ₈; by causing the CPU of the electronic computer to function as a fourth computation unit; and a fifth computation step which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F28] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ³ ^(−1)/r) by causing the CPU of the electronic computer to function as a fifth computation unit.
 24. A non-transitory computer recording medium storing a pairing computation program, wherein an elliptic curve is given as y²=x³+ax, aεF_(p), letting an embedding degree be 8, E be an additive group constituted of rational points on the pairing enabled elliptic curve defined over an extension field F_(p) ⁸, E[r] be a set of rational points having a prime order r, and φ_(p) be a Frobenius endomorphism, a pairing e is defined using G ₁ =E[r]∩Ker(φ_(p)−[1]), G ₂ =E[r]∩Ker(φ_(p) −[p]), as a non-degenerate bilinear map e:G ₁×G₂→F*_(p) ⁸/(F* _(p) ⁸)^(r), the pairing computation program causes an electronic computer which includes a CPU to compute the pairing e(S,Q), by letting SεG₁, QεG₂, χ be a given integer variable, and f_(3χ,S)(Q) be a rational function which is calculated using Miller's algorithm, the order r and a trace t of the Frobenius endomorphism φ_(p) being specified using the integer variable χ as, r(χ)=9_(χ) ⁴+12_(χ) ³+8_(χ) ²+4_(χ)+1, t(χ)=−9_(χ) ³−3_(χ) ²−2_(χ), and a representation of the integer variable χ using p² and p³ with p as a characteristic being p ³≡p²+3_(χ)+1(modr(χ)) the pairing computation program causing the CPU of the electronic computer to function as: an input unit which inputs the integer variable χ, the rational point S, and the rational point Q into respective predetermined registers; a first computation unit which computes f_(3χ,) ^(S) (Q) and f_(3χ,p^3S)(Q); a second computation unit which computes respective rational points p²(S), (3_(χ)+1)S, p²(p³S), (3_(χ)+1)p³S in order using previously obtained results; a third computation unit which respectively computes a value l₅ at a rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (3_(χ)S, S), a value l₆ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p²(S), (3_(χ)+1)S), a value l₇ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (3_(χp) ³S, p³S), and a value l₈ at the rational point Q (x_(Q), y_(Q)) of a straight line passing through rational points (p²(p²S), (3_(χ)+1)p³S); and a fourth computation unit which computes f′_(χ,) ^(S) (Q) using said f_(3χ,) ^(S) (Q) and said f_(3χp^3S)(Q) and said values l₅, l₆, l₇, l₈ as [F29] f′ _(χ,) ^(S) (Q)=(f _(3χ,) ^(S) (Q)·l ₅ ·l ₆)^(p) ³ f _(3χ,p) ₃ _(S)(Q)·l ₇ ·l ₈; and a fifth computation unit which computes the pairing e(S,Q) using said f′_(χ,) ^(S) (Q) as [F30] e(S,Q)=f′ _(χ,) ^(S) (Q)^((p) ³ ^(−1)/r). 